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1.1  SI  OF  SYMBOLS 


The  symbology  for  this  Technical  Note  ls  summarized  in  this  section. 

Definitions  and  acronyms  are  given  in  the  Glossary. 

a  -  Radius  of  the  earth.  The  radius  of  a  sphere  having  the  same 

volume  as  the  <  artii  is  6  3  7  1.22  1  3  km. 

Bn  -  3DNLPI1  box  number. 

Bs  -  SGDB  box  number. 

de  -  Distance  between  grid  points  on  the  earth's  surface. 

d,,'  -  Distance  between  grid  points  on  the  earth's  surface  in  the 

north-south  direction  on  a  Mercator  projection. 

df  -  Distance  between  grid  points  on  the  image  plane. 

dm  -  Distance  between  grid  points  on  the  map. 

d0  -  The  distance  of  181  km.  Among  other  things,  this  distance 

equals  the  whole-mesh  grid  Length  at  60°N  (or  °S)  for 
polar  stereographic  projections,  as  well  as  the  distance 
between  grid  points  on  the  image  plane  for  polar  stereographic 
projections  tine  at  60°N  (or  °S). 

G  -  Grid  designator  for  Mercator  grids. 

H  -  Designator  of  the  hemisphere  on  which  a  polar  stereographic 

projection  is  centered. 

1  -  Indicator  of  position  relative  to  grid  points  in  the  J-x 

direction. 

I'  -  I  index  for  a  second  grid. 

Ipi  -  Total  number  of  points  in  the  east-west  direction  for  a 

Mercator  grid. 

In  -  Number  of  nondup 1 ica t ed  points  in  the  east-west  direction  for 

a  Mercator  grid. 

Ln  -  I  index  for  the  3DNEPH  Box  IJ  convention. 

Inp  -  I  index  of  tiie  North  Pole. 

Ip  -  1  index  of  the  pole. 

Ip'  -  I  index  of  the  pole  for  a  second  grid. 

^s  -  I  index  for  the  SGDB  Box  IJ  convention  as  used  in  this 

Technical  Note. 


V 

^sp 

*w 

V 
12 
4 
18 
:64 

:64  1 
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JE 

JM 
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^2 
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-  I  index  for  the  SGDB  Box  IJ  convention  as  used  in  SYNAPS. 

-  I  index  of  the  South  Pole. 

-  I  index  for  a  whole-mesh  grid. 

-  I  index  for  a  whole-mesh  grid  as  used  in  SYNAPS. 

-  I  index  for  a  half-mesh  grid. 

-  I  index  for  a  quarter-mesh  grid. 

-  I  index  for  an  eighth-mesh  grid. 

-  I  index  for  a  sixty-fourth-mesh  grid. 

-  I  index  for  a  sixty-fourth-mesh  grid  as  used  in  SYNAPS. 

-  Indicator  of  position  relative  to  grid  points  in  the  'y 
d i rec  t ion. 

-  J  index  for  a  second  grid. 

-  J  index  of  the  Equator  for  a  Mercator  grid. 

-  Number  of  points  in  the  north-south  direction  for  a  Mercator 
grid. 

-  J  index  for  the  3DNEPH  Box  IJ  convention. 

-  J  index  of  the  North  Pole. 

-  J  index  of  the  pole. 

-  J  index  of  the  pole  for  a  second  grid. 

-  J  index  for  the  SGDB  Box  IJ  convention  as  used  in  this 
Technical  Note. 

-  J  index  for  the  SGDB  Box  IJ  convention  as  used  in  SYNAPS. 

-  J  index  of  the  South  Pole. 

-  J  index  for  a  whole-mesh  grid. 

-  J  index  for  a  whole-mesh  grid  as  used  in  SYNAPS. 

-  J  index  for  a  half-mesh  grid. 

-  J  index  for  a  quarter-mesh  grid. 

-  J  index  for  an  eighth-mesh  grid. 


-  J  index  for  a  six ty- fourth-mesh  grid  as  used  in  SYNAPS. 

-  Any  integer. 
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-  Mesh  factor. 

-  Mesh  factor  for  a  second  grid. 

-  Map  factor. 

-  Constant  of  the  cone. 

-  Radial  distance  on  the  image  plane  between  the  pole  and  a 
given  latitude  in  multiples  of  grid  length  d^. 

-  Radial  coordinate  of  the  image-plane  polar  coordinate  system 
(  r ,  t) ) . 

-  Scaling  parameter  defined  for  polar  stereographic  projections 
so  that  whole-mesh  grid  points  are  exactly  381  km  apart  at 
60°N  (or  °S)  on  i  lie  earth's  surface,  regardless  of  the 

si aiula rd  latitude. 

-  Wind  component  in  the  ‘x  direction. 

-  Wind  component  in  the  ‘y  direction. 

-  Grid-point  Cartesian  coordinate  in  the  x  direction;  expressed 
in  multiples  of  grid  length  dj  (which  is  a  function  of  mesh 
size)  from  the  y  axis. 

-  X  coordinate  for  a  whole-mesh  grid. 

-  A  coordinate  of  the  image-plane  Cartesian  coordinate  system 
(x,y).  The  distance  x  from  the  y  axis  is  with  respect  to  the 
image  plane,  not  the  surface  of  the  earth. 

-  Grid-point  Cartesian  coordinate  in  the  y  direction;  expressed 
in  multiples  of  grid  length  d£  (which  is  a  function  of  mesh 
size)  from  the  x  axis. 

-  Y  coordinate  for  a  whole-mesh  grid. 

-  A  coordinate  of  the  image-plane  Cartesian  coordinate  system 
(x,y).  The  distance  y  from  the  x  axis  is  with  respect  to  the 
image  plane,  not  the  surface  of  the  earth. 

-  Longitudinal  (east-west)  spacing  between  grid  points. 

-  Latitudinal  (north-south)  spacing  between  grid  points. 

-  Azimuthal  coordinate  of  the  image-plane  polar  coordinate 
system  (  r  ,0  ) . 

-  Longitude. 
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-  Reference  longitude. 

-  Map  scale. 

-  The  constant  3.14159... 

-  Image  scale. 

-  Latitude. 

-  The  standard  latitude  for  a  polar  stereographic  projection. 

-  A  standard  latitude  for  a  Mercator  or  Lambert  conformal 
projection. 

-  A  standard  latitude  for  a  Mercator  or  Lambert  conformal 
projection. 
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SECTION  1.  INTRODUCTION 


The  AFGWC  system  of  maps  and  grids  lias  evolved  over  the  past  20  years. 

The  system  in  its  present  form  is  large  and  complex.  Tarbell  and  Hoke  (1979) 
described  factors  that  have  led  to  the  present  form  of  the  AFGWC 
analysis/ forecast  system.  These  same  factors  have  led  to  the  development  of 
the  system  of  maps  r.nd  grids  used  to  support,  the  AFGWC  analysis/forecast 
system.  These  factors  include: 

a.  AFGWC 1 s  mission  to  provide  worldwide  meteorological  support. 

b.  The  continual  development  of  meteorological  analysis  and  forecast 
models  and  techniques  at  AFGWC  and  other  locations  in  the  meteorological 
community . 

c.  The  specific  hardware  at  AFGWC. 

d.  The  continual  technological  development  of  systems  for  taking 
meteorological  observations,  especially  meteorological  satellites. 

e.  A  greater  availability  of  data  in  the  Northern  Hemisphere  than  in  the 
Southern  Hemisphere. 

f.  A  larger  number  of  operational  requirements  in  the  Northern  Hemisphere. 

g.  The  need  for  grids  compatible  with  those,  of  other  U.S.  numerical 
weather  centers. 

In  view  of  the  above  reasons,  no  master  plan  addressing  the  entire  system 
of  maps  and  grids  was  possible  when  the  system  was  initiated.  Even  today 
changing  operational  requirements  are  reflected  by  the  changing  automated 
analysis/forecast  system.  The  system  of  maps  '•nd  grids  must  change  with  it. 

There  are  two  general  classes  of  environmental  support  provided  at  AFGWC : 
meteorological  and  space  environmental.  Meteorological  support  includes 
observation,  analysis,  and  forecasting  of  the  tropospheric  and  stratospheric 
environment.  Space  environmental  support  includes  observation,  analysis,  and 
forecasting  of  the  space  and  near- space  environment.  Although  some  of  the 
maps  and  grids  are  not  common  to  both,  only  those  used  to  provide 
meteorological  support  will  be  discussed  in  this  Technical  Note. 

The  mathematical  details  of  the  derivations  of  the  various  projections  and 
maps  are  not  presented  here.  Also,  we  do  not  discuss  the  formulation  of 
equations,  such  as  the  equations  of  motion,  on  the  various  projections- 
Excellent  discussions  of  these  topics  may  be  found  in  publications  by  Deetz 
and  Adams  (1945),  Thomas  (1952),  Richardus  and  Adler  (1972),  Gerrity  (1973), 
and  Williamson  (1979),  for  example. 

We  discuss  the  projections  and  maps  most  commonly  used  at  AFGWC  in  Section 
2,  with  the  latitude  and  longitude  conventions  used  in  this  publication 
defined  in  Section  2.1.  Section  3  details  the  AFGWC  grid  systems.  The 
domains  of  the  Northern  and  Southern  Hemispheric  regions  for  identifying  and 
storing  observations  are  depicted  in  Section  4.  Section  5  is  a  glossary  of 
the  most  important  terms  presented  in  this  Technical  Note.  A  list  of 
references  is  presented  in  Section  6. 


SECTION  2.  PROJECTIONS  AND  MAPS 


The  map  has  become  a  basic  tool  of  the  meteorologist.  Saucier  (1955) 
noted  that  meteorological  maps  (charts)  used  to  represent  and  interpret  the 
state  of  the  atmosphere  "must  best  portray  the  space  variations  of  the 
atmospheric  variables  with  consideration  for  convenience  in  plotting  and 
analysis,  for  accuracy  in  representation,  and  for  interpretation".  In  this 
section  we  present  the  mathematical  bases  for  the  polar  stereographic, 
Mercator,  and  Lambert  conformal  projections  and  maps,  which  are  three  popular 
systems  for  satisfying  Saucier's  requirements. 

2.1  CONVENTIONS  FOR  LATITUDE  AND  LONGITUDE 

Before  discussion  of  map  projections,  it  is  necessary  to  specify  the 
conventions  for  locating  points  on  the  earth's  surface  by  latitude  't1  and 
longitude  X.  For  latitude,  Northern  Hemispheric  locations  are  designated  in 
this  Technical  Note  by  °N  increasing  northward  or  by  positive  values 
(without  the  N)  increasing  northward.  Southern  Hemispheric  latitudes  are 
designated  by  °S  increasing  southward  or  by  negative  values  (without  the  S) 
increasing  northwaid.  For  example,  the  South  Pole  is  at  90°S,  or  -90°. 

For  longitude,  Eastern  Hemispheric  locations  are  indicated  by  °E 
increasing  eastward  or  by  positive  values  between  0°  and  180°  (without  the 
E)  increasing  eastward.  Western  Hemispheric  locations  are  indicated  by  °W 
increasing  westward,  negative  values  between  -180°  and  0°  (without  the  W) 
increasing  eastward,  or  positive  values  between  180°  and  360°  (without  the 
W)  increasing  eastward.  As  example,  90°E  is  the  same  as  90°,  and  90°W 
is  the  same  as  -90°  and  270°. 

The  equations  presented  in  subsequent  sections  are  valid  for  those 
conventions  not  using  the  alphabetical  designators  (N,  S,  E,  W).  It  is 
important  to  note  that  in  specific  applications  at  AFGWC ,  latitude-longitude 
conventions  different  from  those  of  this  Technical  Note  might  be  used.  As  an 
example,  in  the  global  Hough  analysis  (HUFANL)  longitude  increases  toward  the 
west. 


2.2  MAP  CONSTRUCTION  AND  SCALE  FACTORS 

As  used  at  AFGWC,  the  term  map  may  be  defined  as  a  two-dimensional 
horizontal  representation  of  the  surface  of  the  earth  with  size  convenient  for 
display,  analysis,  and  interpretation  of  meteorological  data.  The  fundamental 
problem  in  map  construction  is  the  transfer  of  the  earth's  geography  from  its 
actual  three-dimensional  form  to  a  flat  surface  in  such  a  way  as  to  present 
the  earth's  surface  most  advantageously.  Conceptually,  two  transformations 
are  necessary:  first,  the  curved  surface  of  the  earth  must  be  transformed  to 
a  two-dimensional  surface;  second,  this  two-dimensional  surface  must  be 
reduced  to  a  convenient  size. 

The  first  transformation  has  several  steps.  First,  the  irregularly  curved 
surface  of  the  earth  is  transformed  into  a  regularly  shaped  surface,  usually 
the  sphere  or  ellipsoid.  In  subsequent  discussions  in  this  Technical  Note,  we 
will  assume  that  this  first  step  of  the  first  transformation  has  been 
performed  so  that  the  earth  may  be  considered  spherical.  In  the  second  step 


c he  regularly  shaped  three-dimensional  surface  is  transiorined  into  a 
two-dimensional  one.  As  will  be  discussed  in  the  next  section  this  step  can 
he  done  in  three  different  ways.  In  one  way  the  earth's  surface  is  projected 
using  straight  lines,  called  rays,  onto  i  geometric  figure,  such  as  the  cone 
in  Fig.  2.1,  a  cylinder,  or  a  plane.  Then  the  surface  might  be  stretched  to 
produce  some  desirable  quality,  such  as  conservation  of  shape  or  area,  on  the 
resultant  surface  (the  image  surtace).  Finally,  if  necessary,  this  surface  is 
cut  and  unfolded,  a  process  called  ie  ve  1  op  i  ng ,  to  produce  a  planar 
surf  Uvj  o —  l  tie  image  plane  (see  Fig.  2.2).  The  representation  of  the  surface  of 
cue  e.ircti  on  the  image  surface  or  image  plane  is  called  the  pro  gee  t  ion.  The 
t  runs  forma  t  ion  from  curved  eartii  surface  to  image  surface  or  image  plane  is 
described  mathematically,  in  part,  by  tbe  image  scale  o,  which  is  the  ratio  at 
a  point  of  (differential)  distance  on  (lie  image  surface  or  image  plane  to 
(differential)  distance  on  the  surface  ol  the  earth.  The  mathematical 
expression  for  image  scale  varies  according  to  projection  and  is  described  in 
more  detail  in  Sections  2.4.2,  2.5.2,  and  2.6.2. 

The  second  transfonnat ion  consists  of  uniformly  reducing  the  image  plane 
to  a  size  suitable  for  the  purpose  intended  to  produce  a  map ,  or  map  surface. 
This  transformation  is  described  mathematically  by  the  map  scale  p,  which  is 
the  ratio  of  distance  on  the  map  to  distance  on  the  image  surface  or  image 
plane.  The  map  scale  specifies  the  image-plane  reduction  necessary  to  form 
the  particular  map;  for  example,  P  =  1:7,500,000  refers  to  a  map  generated  by 
uniformly  reducing  image-plane  distances  by  a  factor  of  7,500,000.  Maps 
commonly  used  at  AFGWC  vary  in  map  scale  from  1:60,000,000  to  1:3,750,000. 

Both  transformations  can  be  described  bv  one  factor  called  the 
map  factor  m,  where 

(2.1) 

Thus,  in  is  the  ratio  at  a  point  of  (differential)  distance  on  a  map  to 
(differential)  distance  on  the  earth. 


2.3  CLASSIFICATION  OF  PROJECTIONS  AND  MAPS 

As  previously  noted,  a  particular  map  it  selected  so  that  its 
characteristics  best  suit  the  purpose  of  the  user.  Because  the  number  of 
unique  projections  and  maps  is  infinite,  a  general  classification  scheme  is 
needed.  Richardus  and  Adler  (1972)  presented  a  scheme  suggested  by  Goussinsky 
(1951)  in  which  five  general  classes,  which  were  not  mutually  exclusive,  were 
defined.  Each  ulass  was  partitioned  into  mutually  exclusive  subclasses, 
called  varieties.  This  scheme  is  summarized  here  because  it  offers  a  simple, 
yet  comprehensive,  method  of  classifying  projections  and  maps.  The  five 
general  classes  of  the  Goussinsky  scheme  are  nature,  coincidence,  position, 
properties,  and  generation. 

Nature  describes  the  geoine trie  f l gurc  associated  with  the  image  surface. 
The  plane ,  cone ,  and  cylinder  represent  the  only  universally  recognized 
varieties  of  the  nature  class;  all  three  are  associated  with  maps  for 
meteorological  use. 


mm 


Figure  2.2. 
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Transformation  of  image  aur faces  into  image  planes  by 
cutting  and  unfolding  (after  ©Steadier,  1969;  John  Wiley 
and  Sons,  Inc.,  Publishers). 


Figure  2.1. 


Projection  of  the  surface  of  the  earth  onto  a  cone 
(©Strahler,  1969;  John  Wiley  and  Sons,  Inc.,  Publishers) 


Cu  i  nc  idence  refers  t.o  the  L  v ;  -  >  mii.i  In-tween  the  geometric  figure  and 

the  earth.  Three  varieties  have  been  ur  I  •  i .  <  1  (Fig.  2  .  i  )  .  The  tangent  variety 
includes  nlL  projections  in  which  tin  geoine  t  i  i  c  figure  is  tangent  to  the 
earth.  In  the  secant  variety,  the  geu.netri  figure  intersects  the  interior  of 
the  earth.  In  the  po lysupe  r f i c  ia 1  variety,  a  succession  of  different 
geometric  figures  is  used  to  const; a.  !  wv  projection.  Although  the 
polysuperf ic i a  1  variety  is  most  m <  at.  ,  i  he  resulting  maps  are  cumbersome 
and  the  trans  Lormat  ion  equnt  i  oi.s  are  p  it--  complex.  The  tangent  and  secant 
varieties,  on  the  other  hand,  often  pt  oduce  easv-to-use  maps  with  an 
acceptable  amount  of  distortion,  ami ,  consequently,  serve  as  a  basis  for  most 
meteorological  maps. 

Position  re  fers  to  the  alignment  >  the  gf-omet  ric  figure  with  respect  to 
the  axis  of  the  earth  (Fig.  2.4).  In  i  iiorm.il  position,  the  axis  of 
symmetry  of  the  geometric  figure  coin  l.e-a  w,  i  the  rotational  axis  of  the 
earth.  In  the  transverse  position,  th  axis  symmetry  is  perpendicular  to 
the  earth's  axis.  Those  projection-  ,,  winch  Lhe  axis  of  symmetry  is  neither 
perpendicular  to  nor  coincident  with  In  earth's  axis  belong  to  the  oblique 
variety.  The  normal  variety  of  geometric  figure  is  most  widely  used  because  a 
primary  interest  of  the  meteorologist  is  atmospheric  phenomena  that,  in 
general,  move  about  Lhe  earth's  axis. 

There  are  three  mutually  exclusive  vs  i  el  ns  of  Lhe  prope  rty  class. 
liquid  istance  is  defined  as  correct  iepr<  seuini  ion  of  distance  between  two 
points.  This  property  is  limited  to  «  rtai.ii  specified  points  and  cannot  be 
generalized  over  an  entire  image  surl.ne.  In  the  equivalent  variety,  areas  of 
features  on  the  earth  are  preserved  on  :  In-  image  surface,  correct 
representation  of  shapes,  however,  is  sac r s I  iced .  The  conformal  (or 
o rthomo rph ic )  variety  includes  image  surfaces  in  which  the  angles  between 
intersecting  curves  and,  therefore,  t he  shapes  of  small  features  on  the 
earth's  surface  are  preserved.  Also,  a  differentially  small  distance  on  the 
projection  will  correspond  to  a  d  is  la  a  e  earth's  surface  independent  of 

direction  from  the  point.  Con  forma  i  m a  e  used  extensively  in 
meteorology.  Their  accurate  represent af i on  ot  geographical  shapes  facilitates 
the  interpretation  of  meteorological  dun. 

The  gene  rat  ion  class  describes  the  method  used  to  construct  the  image 
surface.  For  the  geometric  (or  £«'  rsjioc  L  i  v<  .)  variety,  the  surface  of  the  earth 
is  simply  projected  using  rays  onto  a  geometric  figure.  The  conventiona 1 
variety  includes  projections  generated  puioly  mathematically  without  regard  to 
projection  using  rays.  The  seim-gcomet  ri c  variety  is  a  combination  of  the 
geometric  and  conventional  varieties. 

The  Gouss  inslcy  classification  scheme  is  summarized  in  Table  2.1.  As  is 
readily  apparent,  there  is  a  large  number  of  unique  possibilities. 

Projections  used  as  a  basis  for  meteorological  maps  represent  only  a  small 
subset  of  this  number.  Ac  AFG'WC ,  extensive  use  is  made  of  polar  stereographic 
and  Mercator  maps,  produced  from  polar  stereographic  and  Mercator  projections, 
respectively.  In  addition,  some  specialized  support  is  provided  using  Lambert 
conformal  maps.  The  Goussinsky  classification  as  it  relates  to  AFGWC 
projections  and  maps  is  shown  in  Table  2.2 .  The  three  projections  and  maps 
most  commonly  used  in  meteorological  applications  will  be  described  in  detail 
in  thi'  remainder  of  this  chapter. 
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Table  2.1.  Summary  of  the  Goussinsky  Classification  Scheme  (Goussinsky,  1951). 


CLASS 

VARIETIES 

Nature 
Coinc idence 
Position 
Property 
Gene  ration 


Planar,  Conical,  Cylindrical 
Tangent,  Secant,  Polysuperf ic ia 1 
Normal,  Transverse,  Oblique 
Equidistant,  Equivalent,  Conformal 
Geometric,  Conventional,  Semi-geometric 


Table  2.2.  Goussinsky  Classification  Scheme  for  AFGWC  Projections  and  Maps. 


PROJECTION/MAP 

NATURE 

COINCI¬ 

DENCE 

POSITION 

PROPERTY 

GENERATION 

Northern  Hemispheric 
Polar  Stereographic 

Planar 

Secant 

Normal 

Conf orma 1 

Geometric 

Southern  Hemispheric 
Polar  Stereographic 

Planar 

Secant 

Normal 

Conf orma 1 

Geometric 

Mercator 

Cylin¬ 

drical 

Secant 

Normal 

Conf orma 1 

Semi- 

geometric 

Lambert  Conformal 

Conical 

Secant 

Normal 

Conformal 

Conventional 

2.4  POLAR  STEREOGRAPH! C  PROJECTIONS  AND  MAPS 


The  secant  polar  stereographic  projection  is  used  extensively  at  AFGWC . 
This  projection  is  generated  geometrically  (that  is,  using  rays)  by 
positioning  a  secant  plane  normal  to  the  earth's  axis  at  a  standard  latitude. 
(Fig.  2.5a  illustrates  the  geometry  for  the  Northern  Hemispheric  polar 
stereographic  projection.)  The  resulting  projection  has  the  conformal 
property.  In  the  special  case  that  the  standard  latitude  is  a  pole  (4>f  = 

90°N  and  90°S  for  Northern  Hemispheric  and  Southern  Hemispheric  polar 
s te reographic  projections,  respectively),  a  tangent  polar  stereographic 
projection  is  formed  (Fig  2.5b).  The  tangent  projection  is  used  at  some 
meteorological  facilities  although  AFGWC  is  not  one  of  them. 

In  the  following  discussion,  the  equations  for  the  secant  polar 
stereographic  projection  and  map  will  he  presented.  The  equations  for  the 
tangent  case  are  obtained  using  $r  =  90°N  or  °S . 

2.4.1  GENERATION  OF  THE  SECANT  POLAR  STEREOGRAPHIC  MAP 

The  secant  polar  stereographic  map  is  formed  in  the  following  manner. 
First,  a  plane,  which  is  the  geometric  figure  for  this  projection,  is  passed 
through  the  earth  (perpendicular  to  the  earth's  axis)  at  a  standard  latitude, 
or  true  latitude ,  4>0  >  as  shown  in  Fig.  2.5a.  At  AFGWC,  4^  =  60°N  is  used 
for  the  Northern  Hemispheric  polar  stereographic  projection  and  4*  =  60°S 
for  the  Southern  Hemispheric  polar  stereographic  projection.  Next,  for  each 
latitude  4>  a  straight  line,  or  ray,  is  drawn  through  the  pole  (point  P)  of  the 
opposite  hemisphere  and  the  surface  of  the  earth  (point  E).  This  ray  (ray  PIE 
in  Fig.  2.5a)  projects  the  surface  of  the  earth  onto  the  plane,  the  image 
plane  and  image  surface  in  this  case,  at  point  I.  After  the  desired  points  on 
the  earth's  surface  are  projected  the  image  plane  is  reduced  to  a  map  of 
manageable  size  by  means  of  the  map  scale  I*  (see  Section  2.2). 

The  image  plane  (image  surface)  and  map  (Fig.  2.6)  are  conformal  and  have 
the  following  characteristics: 

a.  Lines  of  constant  latitude  are  concentric  circles  centered  at  the 

po  le . 

b.  Lines  of  constant  longitude  (meridians)  are  depicted  as  straight 
lines  radiating  from  the  pole.  These  radials  are  spaced  at  the  same  angular 
increment  as  meridians  on  the  earth. 


2.4.2  POLAR  STEREOGRAPHIC  IMAGE  SCALE 

As  noted  previously,  image  scale  accounts  for  the  transformation  from  the 
earth's  surface  to  the  image  surface  or  image  plane  and  is  the  ratio  of 
distance  on  the  image  surface  or  image  plane  to  distance  on  the  earth's 
surface.  For  a  polar  stereographic  projection  true  at  latitude  <t>Q,  the  image 
scale  o  is  given  by 

1  +  H  sin  if 

nU)  =  — - -  ,  (2.2) 

1  +  H  s  i  n  <j' 

where  <{>  is  latitude.  Eq.  (2.2)  is  valid  for  both  Northern  and  Southern 
Hemispheric  polar  stereographic  projections,  with 
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+  ]  ,  Northern  Hemispheric  project,  ion, 
-1,  Southern  liemi :  pheri  c  pro,]'  ct,  i  on  . 


(2.3) 


Note  that  Southern  Hemispheric  locations  (<f)  <  0°)  may  exist  on  a  Northern 

Hemispheric  projection  (H  =  +1),  and  vice  versa.  The  latitudinal  variation  of 

image  scale  is  shown  in  Fig.  2.7.  The  solid  curve  is  the  secant  case 

for  4>o  =  60°N  or  °S  and  the  dashed  curve  is  the  tangent  case 

((f)  =  90°N  or  °S).  Earth  distance  is  preserved  on  the  image  plane  (image 

surface)  only  at  the  true  latitude  tf>  ,  where  o  =  1.  Poleward  of  latitude  <|> 

.  .  .  u  o 

the  distance  between  two  points  is  larger  on  the  earth  than  on  the  image  plane 

(image  surface),  whereas  equatorward  of  <!>  the  distance  is  smaller  on  the 

earth. 


2.4.3  IMAGE-PLANE  COORDINATE  SYSTEMS 


An  image-plane  Cartesian  coordinate  system  (x,y)  can  be  defined  for  the 
polar  stereographic  image  planes.  The  orientation  of  the  coordinate  Sxes  for 
the  Northern  and  Southern  Hemispheres  is  shown  in  Figs.  3.6  and  3.7, 
respectively.  In  both  hemispheres  the  origin  is  the  pole,  the  positive  x  axis 
coincides  with  the  10°E  meridian,  and  the  positive  y  axis  coincides  with  the 
100°E  meridian.  These  orientations  result  in  a  right-handed  system  for  the 
Northern  Hemispheric  projection  and  a  left-handed  system  for  the  Southern 
Hemispheric  projection.  The  image-plane  Cartesian  coordinates  of  any  point  on 
the  earth's  surface  are  given  by 


x  =  o  n  cos  I  cos(A-A^) 


(2.4) 


y  =  o  a  cot; 


t;  (f  s>in(A-A  ) 


(2.3) 


where  o  is  the  image  scale,  ij>  is  the  latitude,  A  is  the  longitude,  A  is  the 
reference  longitude  of  the  positive  x  axis  (at  AFGWC ,  A  =  10°E),  an§  a  is 
the  earth's  radius.  In  this  Technical  Note  we  use  the  radius  of  a  sphere 
having  the  same  volume  as  the  earth  for  a, 


a  =  6371.2213  km 


(2.6) 


Certain  numerical  models  at  AFGWC  may  use  slightly  different  values  for  a. 

A  corresponding  image-plane  polar  coordinate  system  (r,0)  has  also  been 
defined.  As  shown  in  Figs.  3.6  and  3.7,  radial  distance  r  increases  with 
distance  from  the  pole.  Also,  azimuth  0  increases  from  the  reference 
longitude,  \  =  10°E,  in  the  counter-clockwise  direction  on  the  Northern 

Hemispheric  projection  and  the  clockwise  direction  on  the  Southern  Hemispheric 
projection.  The  polar  coordinates  on  the  image  planes  are  given  by 


r  =  o  it  cos  .(> 


(2.7) 


0  =  \  -  A 


(2.8) 


.  i; 
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Figure  2.7. 


The  latitudinal  variation  of  image  scale  o  for  polar 
stereographic  projections  true  at  60°N  or  °S  (solid 
line)  and  90°N  or  °S  (dashed  line).  Not  plotted  is  the 
image  scale  for  these  projections  extending  into  the 
opposite  hemisphere. 


2.5  MERCATOR  PROJECTIONS  AND  MAPS 


A  secant  Mercator  projection  is  used  at  AFGWC  for  certain  tropical 
applications.  Although  in  our  discussion  we  employ  the  cylinder  as  the 
geometric  figure  used  in  the  construction  of  the  Mercator  projection,  this 
projection  was  originally  derived  without  consideration  of  geometric  figure  by 
Gerhardus  Mercator  for  its  navigational  benefits  (Deetz  and  Adams,  1945). 

This  conformal  projection  can  be  generated  semi-geomet rica 1 ly  by  projecting 
the  earth's  surface  onto  a  normally  positioned  cylinder  that  intersects  the 
earth  at  standard  latitudes  ^  and  1»2  (Fig.  2.8a)  and  then  shrinking  t he 
result  in  the  north-south  direction.  For  the  case  where  'h  =  d*2  =  0° 

(Fig.  2.8b)  the  tangent  Mercator  projection  is  obtained. 

In  the  following  discussion,  the  equations  of  the  secant  Mercator 
projection  and  map  will  be  presented.  Equations  for  the  tangent  case  are 
obtained  using  !]  =  ^  =  0°. 

2.5.1  GENERATION  OF  THE  SECANT  MERCATOR  MAP 

The  secant  Mercator  map  can  be  derived  in  the  following  manner.  First,  a 
cylinder,  which  is  the  geometric  figure  for  this  projection,  is  aligned  so 
that  its  axis  coincides  with  the  earth's  axis  and  intersects  the  earth's 
surface  at  latitudes  ij  and  1*2 .  For  AFGWC  projections,  ‘J'p  =  22.5°N 
and  4>2  =  22.5°S.  Secondly,  for  each  latitude  'l1  a  ray  is  drawn  through  the 
center  of  the  earth  (point  0)  and  the  earth's  surface  (point  E).  This  ray 
(ray  OEI  in  Fig.  2.8a)  projects  the  surface  of  the  earth  onto  the  geometric 
figure,  the  cylinder,  at  point  I.  In  Fig.  2.8a,  ray  01 ' E '  projects  the 
Equator  onto  the  cylinder.  After  the  desired  points  on  the  earth's  surface 
are  projected,  the  resultant  is  shrunk  in  the  north-south  direction  to  obtain 
an  image  surface  that  is  conformal.  Next,  the  image  surface  is  developed  into 
the  image  plane  by  cutting  the  image  surface  along  a  reference  meridian  and 
unrolling.  Finally,  the  image  plane  is  reduced  to  a  map  of  manageable  size  by 
means  of  the  map  scale  •* . 

On  the  image  plane  and  map  (Fig.  2.9),  latitude  circles  are  depicted  as 
parallel  straight  lines  spaced  as  a  function  of  latitude.  Meridians  are 
depicted  as  equally  spaced,  parallel  straight  lines  and  are  oriented 
perpendicular  to  the  latitude  lines.  Spacing  between  meridians  on  the  image 
plane  is  equal  to  spacing  between  meridians  on  the  earth  at  the  standard 
latitudes.  Loxodromes  (lines  of  constant  bearing)  are  straight  lines  on  the 
image  plane  and  map. 

2.5.2  MERCATOR  IMAGE  SCALE 

For  Mercator  projections  true  at  latitudes  'J' j_  and  Ty ,  image  scale  is  a 
function  of  latitude  >(i  and  is  given  by 

o(c|.)  =  cor.  I  nee  .|>  =  cos  f  sec  <|>  .  (2.9) 

The  latitudinal  variation  of  image  scale  is  shown  in  Fig.  2.10  where  the  solid 
curve  represents  the  secant  projection  (for  =  22.5°N,  <|> ^  =  22.5°S) 
and  the  dashed  curve  represents  a  tangent  projection  (<|<|  =  !•,  =  0°).  In 
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tliLs  secant  case,  earth  distance  ia  preserved  on  the  image  plane  (image 
surface)  at  22.5°N  and  22.!3°S  (where  0  =  1).  Equatorward  of  these 
latitudes,  o  is  less  than  ur,  ;  so  that  the  distance  between  two  points  is 
greater  on  the  earth  than  on  the  image  plane  (image  surface) ,  whereas  poleward 
of  4' i  and  4>2  the  distance  is  smaller  on  the  earth. 

2.5.3  IMAGE-PLANE  COORDINATE  SYSTEM 


A  right-handed  image-plane  Cartesian  coordinate  system  (x,y)  is  defined 
for  the  Mercator  projection.  Orientation  of  the  coordinate  axes  is  shown  in 
Fig.  3.24.  The  (x,y)  coordinates  of  any  point  on  the  image  plane  are  given  in 
terms  of  the  poinc's  latitude  and  longitude  by  the  following  equations: 


I  Mi)1 


(  \-A  }  con  J 
<’  I 


(\-\  )  con  .f. 


a  and  <  in  decree:;, 

o 


\  ani  a  in  radians, 

o 


(2.10) 


y  =  a  ror.  <!'  i’n[  fan  'i(  +  ) ) 

1 


a  cos  )> 

l 


in 


1+r.  1  mj> 

COS(J> 


(2.11) 


In  these  equations  A  is  longitude  and  Aq  is  the  reference  longitude.  At 
AFGWC ,  is  0°  for  Mercator  projections.  Also,  4  is  latitude,  cos  4^  (  = 
cos  4' )  ls  the  cosine  of  the  true  latitude,  11  =  3.14159...,  and  £  is  the 
radius  of  the  earth.  By  eqs.  (2.10)  and  (2.11),  x  increases  to  the  east  and  y 
increases  to  the  north. 


2.6  LAMBERT  CONFORMAL  PROJECTIONS  AND  MAPS 

Lambert  conformal  projections  of  the  secant  variety  are  occasionally  used 
at  AFGWC  for  areas  centered  in  the  middle  latitudes,  especially  when  the  area 
has  a  longer  east-west  extent  than  north-south.  To  aid  in  the  general 
understanding  of  the  characteristics  of  the  Lambert  conformal  projection,  in 
the  discussion  to  follow  we  project  the  earth's  surface  onto  a  normally 
positioned  cone  that  intersects  the  earth  at  standard  latitudes  4^  and 
'(>2  (Fig.  2.11a).  Then  we  shrink  or  stretch  the  resultant  in  both  the 
north-south  and  east-west  directions.  In  actuality,  however,  the  secant 
Lambert  conformal  projection  is  derived  mathematically  (the  conventional 
generation  of  the  Goussinsky  classification  scheme).  A  Lambert  conformal 
projection  of  the  tangent  variety  (true  at  4i  =  't^)  on  the  other  hand,  can 
be  constructed  semi-geome t rica 1 ly . 

In  the  following  discussion,  equations  for  the  secant  projection  and  map 
will  be  presented.  Equations  for  the  tangent  case  are  obtained  with 
4j  =  42,  except  as  noted  in  Section  2.6.2. 
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2.6.1  GENERATION  OF  THE  SECANT  LAMBERT  CONFORMAL  MAP 


The  secant  Lambert  conformal  map  can  be  formed  in  the  following  manner. 
First,  a  cone,  which  is  the  geometric  figure  for  this  projection,  is 
positioned  so  that  its  apex  is  directly  over  one  of  the  earth's  poles  and  its 
sides  pass  through  standard  latitudes  <!>l  and  $2  (Fig.  2.11a).  Next,  for 
each  latitude  4*  the  surface  of  the  earth  (point  E)  is  projected  onto  the  cone 
(the  geometric  figure)  at  point  1  by  drawing  a  ray  from  the  center  of  the 
earth  (point  0)  (Fig.  2.11a).  After  the  desired  points  on  the  earth's  surface 
are  projected,  the  resultant  is  shrunk  or  stretched  in  both  the  north-south 
and  east-west  directions  to  obtain  an  image  surface  that  is  conformal.  Then 
the  image  surface  is  developed  into  the  two-dimensional  image  plane  by  cutting 
the  image  surface  along  a  meridian  and  unrolling  it.  Finally,  the  image  plane 
is  reduced  to  a  map  of  manageable  size  by  means  of  the  map  scale  h . 

The  image  plane  and  map  (Fig.  2.12)  are  conformal  everywhere  except  at  the 
pole.  Meridians  are  depicted  as  straight  lines  converging  at  the  pole. 
Latitude  circles  appear  as  arcs  of  concentric  circles  centered  at  the  pole. 

In  addition,  straight  lines  on  the  image  plane  and  map  approximate  great 
circles.  For  many  meteorological  applications,  only  that  part  of  the 
projection  in  the  region  of  interest  is  displayed. 


2.6.2  LAMBERT  CONFORMAL  IMAGE  SCALE 

For  secant  Lambert  conformal  projections  true  at  ^  and  (J)2 ,  image 
scale  o  is  a  function  of  latitude  and  is  given  by 


where 


:on  >|>  tnnli(90Ol,-<J'  ) 
con  4’ 


»U)  =- 


roc  4, 

1  -n 

1  + 

H  t-.i  n  4.  ^ 

.•-■on  4 

1  + 

11  r in  4 
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11  sin  4> 

con  <|- 
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H  is  gLven  by  eq.  (2.3)  and  n  is  called  the  constant  of  the  cone.  Note  that  n 

is  a  function  of  <f>  ^  and  <f>2 ,  the  latitudes  at  which  the  cone  intersects  the 

earth.  For  <(>]_  =  60°N  (or  °S)  and  <j>2  =■  30°N  (or  °S),  n  =  0.7156. 

In  the  tangent  case,  where  4>i  =  <l'2 1  eq.  (2.13)  is  indeterminate.  In  this 
case , 

n  =  (1  sin  <|’^  =  H  sin  f>2  .  (2.14) 

For  a  cone  tangent  to  the  earth  at  45°N  (or  °S)  ,  n  =  0.7071.  Latitudinal 
variation  of  Lambert  conformal  image  scale  is  shown  in  Fig.  2.13.  The  solid  curve 
represents  a  secant  case  in  which  =  60°N  (or  °S)  and  =  30°N  (or  °S); 

the  dashed  curve  represents  a  tangent  case  in  which  i))j  =  4>2  =  45°N  (or  °S). 

The  small  variation  of  a  in  the  middle  latitudes  is  a  major  reason  why  the  Lambert 

conformal  projection  is  used  in  mid-latitudinal  applications. 

2.6.3  IMAGE-PLANE  COORDINATE  SYSTEMS 

An  image-plane  Cartesian  coordinate  system  (x,y)  is  defined  tor  the  Lambert 
conformal  projections.  Orientation  of  the  coordinate  axes  is  shown  in  Fig.  2.12 
for  a  Northern  Hemispheric  example.  For  both  Northern  and  Southern  Hemispheric 
projections  the  (x,y)  coordinates  of  any  point  on  the  image  plane  are  given  by 


x  =  r  cos  0 


y  =  r  sinfi  , 

where 

o  a  cos  0 

r  =  - - - 

0  =  n  (  A  -  A  q  )  . 

Here  (r,0)  are  the  polar  coordinates  on  the  image  plane  with  r  being  the  radial 
distance  from  the  pole  and  9  being  the  azimuth  (Fig.  2.12).  Also,  A  is 
longitude,  Aq  is  the  reference  longitude  and  is  offset  180°  longitude  from  the 
meridian  of  the  cut  used  in  forming  the  image  plane  from  the  image  surface,  tj)  is 
latitude,  a_  is  the  radius  of  the  earth,  and  a  and  n  are  given  by  eqs.  (2.12)  and 
(2.13),  respectively. 

For  the  special  case  when  n  =  1,  the  polar  stereographic  projection  results. 
The  Mercator  projection  is  the  special  case  when  (ji  = 
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SECTION  3.  GRIDS 


The  advent  of  high-speed  digital  computers  in  the  1940s  offered  new 
opportunities  in  meteorology.  Fot  the  first  time  it  became  feasible  on  a 
near-rea 1-t ime  basis  to  produce  numerical  weather  forecasts  based  on  the 
hydrodynamic  equations.  In  addition,  miputers  could  be  used  to  display  and 
analyze  large  quantities  of  data,  thereby  reducing  the  amount  of  time  and 
effort  required  by  the  meteorologist  to  assimilate  meteorological  information. 


To  facilitate  manipulation  of  data  on  these  computers,  grid  systems  were 
developed.  The  term  gr id  may  be  defined  as  an  array  of  points.  For 
me t eo ro log i ca 1  applications  these  points  are  located  at  the  intersections  of 
uniformly  spaced  parallel  and  perpendicular  lines  (Fig.  3.1).  In  meteorology 
(as  well  as  other  fields  such  as  geography  and  navigation),  this  array  is 
given  meaning  when  it  is  superimposed  on  a  map  (Fig.  3.2)  so  that  individual 
grid  points  represent  points  on  the  surface  of  the  earth.  An  indexing 
convention  called  the  (1,J)  indexing  convention  has  been  developed  to  identify 
individual  grid  positions.  An  example  of  the  (1,J)  indexing  convention  used 
for  most  grids  at  AFGWC  is  shown  in  Fig.  3.3. 


Of  course,  observations  of  the  atmosphere  normally  are  not  taken  at  grid 
points,  but  instead  at  irregularly  spaced  points.  The  process  whereby  the 
data  are  interpolated  from  the  observational  points  to  the  grid  points  is 
called  analysis.  At  AFGWC  analysis  of  meteorological  data  is  performed  by 
numerous  models — the  Hough  global  spectral  analysis  model  (HUFANL)  and  the 
Three-dimensional  Nephanalysis  Model  (3DNEPH)  are  important  examples.  At 
AFGWC,  world-wide  observations  from  Land  surface  stations  and  radiosondes,  for 
example,  are  stored  individually  in  data  bases  separate  from  the  gridded 
data.  (Appendix  A  provides  the  numbered  geographic  regions  in  which  the 
observations  are  located  globally  by  the  World  Meteorological  Organization 
(WMO)  and  are  stored  in  the  AFGWC  data  base.) 


At  AFGWC,  reference  grids  have  been  defined  to  provide  a  basis  for  the 
many  grids  used.  For  a  given  image  surface  (such  as  a  Northern  Hemispheric  or 
Southern  Hemispheric  polar  stereographic  projection),  a  common  geographical 
region  is  specified  for  all  reference  grids.  The  shortest  distance  between 
adjacent  grid  points  on  the  polar  stereographic  Whole-mesh  Reference  Grids  is 
referred  to  as  the  whole-mesh  grid  length.  All  f iner-resolution  grids 
superimposed  on  a  given  image  surface  are  defined  relative  to  the  Whole-mesh 
Reference  Grids.  For  example,  distance  between  half-mesh  grid  points  is 
one-half  the  distance  between  whole-mesh  grid  points;  distance  between 
quarter-mesh  grid  points  is  one-fourth  that  between  whole-mesh  grid  points. 

The  relationship  between  a  whole-mesh  grid  and  its  half-n.esh  and  quarter-mesh 
counterparts  is  shown  in  Fig.  3.4. 


A  number  of  smal  ler-area  grids  that  are  subsets  of  either  the  Whole-mesh 
Reference  Grids  or  a  finer-mesh  version  of  these  reference  grids  are  also 
defined.  These  areal  subsets,  which  cover  only  a  specific  area  of  interest, 
are  used  to  economize  computer  resources. 


In  this  section,  AFGWC  grid  systems  will  be  described.  Included  are  polar 
stereographic,  Mercator,  and  latitude-longitude  grid  systems.  For  each  of 
these  three  systems  a  description  is  given  for  the  reference  grid,  finer-mesh 
versions,  and  areal  subsets,  when  applicable.  This  information  has  been 
discussed,  in  part,  in  publications  by  Headquarters  Third  Weather  Wing  (1962), 
Johnson  (1977),  and  Headquarters  AFGWC  (1978). 


Figure  3.3. 


An  example  of  the  (I,J)  indexing  convention  used  for  most 
grids  at  AFGWC .  In  this  example,  the  column  index  I 
increases  from  left  to  right,  whereas  the  row  index  J 
increases  from  top  to  bottom. 


QUARTER-MESH  POINT 
HALF-MESH  POINT 

WHOUl-MESH  point 


Relationship  between  whole-mesh,  half-mesh,  and 
quarter-mesh  grid  points.  Note  that  each  whole-mesh  point 
is  also  a  half-mesh  and  quarter-mesh  point.  Similarly, 
each  half-mesh  point  is  also  a  quarter-mesh  point. 
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It  is  important  to  note  on  the  following  figures  that  lines  establishing 
boundaries  of  the  various  grids  pass  through  the  grid  points,  not  between  them. 

3.1  POLAR  STEREOGRAPHIC  GRIDS 

Polar  stereographic  grids  are  used  extensively  at  AFGWC  for  numerical 
computations  and  automated  displays.  As  the  name  implies,  these  grids  are 
based  on  polar  stereographic  projections.  Two  reference  domains  have  been 
defined  at  AFGWC — one  centered  on  the  Northern  Hemisphere  and  the  other 
centered  on  the  Southern  Hemisphere.  The  domains  are  based  on  secant  polar 
stereographic  projections  true  at  60°N  and  °S,  respectively.  (The  secant 
polar  stereographic  projection  is  discussed  in  Section  2.4).  For  both  the 
Northern  and  Southern  hemispheric  areas,  reference  grids  of  various 
resolutions  have  been  defined.  Polar  stereographic  grids  for  numerical 
analysis  and  forecasting  at  AFGWC  are  subsets  of  these  reference  grids. 


3.1.1  GENERALIZED  EQUATIONS 


For  polar  stereographic  projections  we  define  a  whole-mesh  grid  as  one  in 
which  the  grid  points  are  exactly  381  km  apart  at  60°N  or  °S  on  the 
earth's  surface.  Such  a  grid  could  be  constructed,  for  example,  by 
superimposing  an  array  of  points  one-half  inch  apart  on  a  polar  stereographic 
map  with  map  scale  1:30,000,000  and  true  at  60°N  (or  °S).  In  this 
Technical  Note  we  represent  the  distance  381  1cm  by  dQ;  that  is, 

d  =  381  km  .  (3.D 

o 

The  distance  between  whole-mesh  grid  points  on  the  earth's  surface  is  a 
function  of  latitude  and  is  named  the  whole-mesh  grid  length  Tor  grid 
interval ,  grid  increment,  grid  unit,  etc.).  An  increasingly  popular  name  for 
this  unit  is  the  bed ient  (Stackpole,  1978),  which  is  named  for  H.  A.  Bedient, 
who  was  instrumental  in  developing  the  map  and  grid  systems  of  the  Joint 
Numerical  Weather  Prediction  Unit  in  the  1930s.  Other  grid  lengths  are  simply 
fractions  of  the  whole-mesh  grid  length — half  mesh,  quarter  mesh,  eighth  mesh, 
and  sixty-fourth  mesh,  for  examples. 


The  distance  d^  between  grid  points  on  the  image  surface  is  given  by 

d.  =  d  0  M  ,  (3.2) 

1  u 

where  M  is  the  mesh  factor,  which  is  defined  as  the  ratio  of  the  grid  length 
to  the  whole-mesh  grid  length  on  the  image  plane.  (Therefore,  for  a 
whole-mesh  grid,  M  =  1,  and  for  a  half-mesh  grid,  M  =  1/2.)  The  scaling 
parameter  S,  where 


I  +  in  60° 


(3.3) 


is  defined  so  that  whole-mesh  grid  points  are  exactly  381  kin  apart  at  60°N 
(or  °S)  on  the  earth's  surface,  regardless  of  the  standard  latitude  <|>0.  For 
applications  ac  AFGWC,  (j>  is  60°N  (and  °S)  so  that  S  =  1.  H  designates 
the  hemisphere  of  the  projection  (see  eq.  (2.3))  and  dG  equals  381  km. 
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By  definition  of  image  scale  o  (Section  2.2)  and  eqs.  (2.2)  and  (3.2),  the 
distance  de  between  grid  points  on  the  earth's  surface  is 

d  C  M 

-  ,  (3.4) 

n 

whe  re 

1  +  !i  r.in* 

nU)  = - ,  (3.5) 

1  +  II  r.  i  ti  i(> 

which  is  valid  for  both  Northern  and  Southern  Hemispheric  projections. 

By  substitution  of  eqs.  (3.3)  and  (3.5)  into  eq.  (3.4),  dg  can  be  shown  to 
be  independent  of  standard  latitude.  For  a  whole-mesh  grid  (M  =  1)  on  a 
Northern  Hemispheric  map  (H  =  +1),  dg  is  381  km  at  60°M ,  as  expected.  A 
plot  of  dg  aa  a  function  of  latitude  is  shown  in  Fig.  3.5  for  a  whole-mesh 
grid. 

By  the  definition  of  map  scale  |i  (Section  2.2)  and  eq.  (3.2),  the  distance 
dm  between  grid  points  on  the  map  is  given  by 

ti  =  |i  d.  =  u  d  f>  M  .  (3.6) 

m  i  o 

Thus,  for  a  whole-mesh  (M  =  1)  polar  stereographic  grid  with  map  scale  P  = 
1:15,000,000  and  true  at  60°N  or  °S,  dm  =  2.54  cm  =  1  inch. 

The  (I,J)  indexing  convention  for  the  grids  identifies  positions  relative 
to  grid  points.  The  convention  is  left-handed  for  the  polar  stereographic 
grids:  the  column  index  I  increases  from  1  on  the  left  of  the  grid  toward  the 
right,  whereas  the  row  index  J  increases  from  1  at  the  top  of  the  grid  toward 
the  bottom.  Figs.  3.6  and  3.7  provide  examples  for  the  Northern  Hemispheric 
and  Southern  Hemispheric  Whole-mesh  Reference  Grids,  respectively. 

The  grid-point  Cartesian  coordinates  (X,Y)  for  a  point  are  expressed  in 
multiples  of  grid  length  d^  (which  is  a  function  of  mesh  factor  M)  from  the 
pole  according  to 


X  =  I  -  I 

p 

•> 

(3.7) 

Y  =  -H  (J  -  J  ) 

•> 

(3.8) 

where  Ip  and  Jp  are  the  I  and  J  coordinates,  respectively,  of  the  pole  on 
which  tne  projection  is  centered.  The  point  (lp,Jp)  may  actually  be 
outside  the  domain  of  the  grid,  as  with  the  U.S.  Boundary-layer  Model  grid, 
tor  example.  H  is  + 1  for  Northern  Hemispheric  grids  and  -1  fcr  Southern 
Hemispheric  grids  (eq.  2.3).  Eqs.  (3.7)  and  (3.8)  are  valid  for  any  point  on 
the  earth.  Examples  for  the  Northern  Hemispheric  and  Southern  Hemispheric 
Whole-mesh  Reference  Grids  are  given  in  Figs.  3.6  and  3.7,  respectively. 
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GRID  LENGTH 
KILOMETERS 


Figure  3.5. 


LATITUDE  (  ^  w  °S) 


The  grid  length  de  on  the  earth's  surface  as  a  function 
of  latitude  4>  for  a  whole-mesh  grid  on  a  Northern 
Hemispheric  or  Southern  Hemispheric  polar  stereographic 
projection. 
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The  Northern  Hemispheric  Whole-mesh  Reference  Grid.  The 
domain  for  this  polar  stereographic  projection  extends  into 
the  Southern  Hemisphere,  which  is  not  plotted.  Also 
included  are  the  axes  for  the  (l,J)  indexing  convention, 
the  image-plane  Cartesian  coordinates  (x,y),  the  grid-point 
Cartesian  coordinates  (X,Y),  and  the  image-plane  polar 
coordinates  (r, 0).  The  w  subscripts  indicate  values  for 
the  whole-mesh  grid.  The  longitude  is  A  and  AQ  is  the 
reference  longitude  ( 10°E  at  AFGWC ) .  The  distances  for 
the  ( x , y 1  coordinates  are,  by  definition,  with  respect  to 
the  image  plane.  The  whole-mesh  grid  spacing  is  displayed 
in  the  upper-left  corner. 
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i  ’-f) 


(2  4  (  Cartesian  coordinates  (*,y),  discussed  in  Section  2.4.3  (eqs, 

grid-poxnt  car  es?  ?  ^  "  "T*  3’6  and  3'7>  «n  be  calculateo  from  ’ 

d.  the  d  st  !  coordinates  (X,Y)  by  multiplying  each  (X,Y)  value  by 

That  is  Ce  betv;een  8rid  points  on  the  image  plane  (see  eq.  (3.2)). 


x  =  d . 


(3.9) 


•v  "  <Ji  Y  ■  (3.10) 

From  eqs.  (2 . 2 )- (2 . 5 )  ,  (3.2)  (3  9i  n  mi  ,  „  •  .  . 

the  latitude  of  anv  end  '  t  \  ,10)’  and  trigonometric  identities 

y  grid  point  may  be  calculated  as  follows: 


o  re::  i  rH 


1  +  H  r,  i  n  tp 


d  SM J  1  +  H  +  X?  +  Y;' 

o  1  K  UJ 


-  X'  +  V 


.  (3.11) 


v:™  "d°s at . 

^^rrrr/s  *«•- in  «•*>.- Ms  Wi-v 

respectively.  T&  a^c  sin"  WtL"  ^s'prin^^e  ’  vil^es  ’  frot  -^o  +90o. 
can  be™  ilcllsilT °‘  ^  ^  POi"' 


A0  +  arccos  rr  Tri  a  Y  >  0 
X  +  Y 


^  ~  A  -  n  rcco:; 


X  +  Y 


?)‘2  ,  Y  <  0 


(3.12) 


(by  definition) 


Y  =  0 


Once  again,  X  and  Y  are  defined  in  eqs.  (3  7)  and  <  i  Si 

reference  longitude,  is  10°E  at  AFGWC  ThI  •  '  .*  N?te  that  \Q,  the 

values  from  0°  to  180°.  arccosine  function  has  principle 


The  (I,J)  coordinates  of  any  point  may  be  expressed  in  terms  of  the 
point's  latitude  and  longitude  according  to 


1  -  r„  +  t-o:i(  .\-.\ 


J  -  <1  -  H  K  r.  i  n  ( A  -  a  ) 

f)  M  o 


(3.13) 

(3.14) 


where 


ti  _  n.  u  cor.  ‘|>  n  o  cor  i|> 
HM  ~ 


.1  M 
o 


(3. IS) 


R(^  is  the  radial  distance  on  the  image  plane  between  the  pole  and  latitude  <1> 
in  multiples  of  grid  length  d^  (see  eq.  (3.2)).  The  radius  of  the  earth  is 
represented  by  ji,  the  mesh  factor  is  M,  longitude  is  A,  reference  longitude 
is  A0)  the  p  subscripts  designate  the  coordinates  of  the  pole,  dQ  is  381 
km,  and  a  is  given  by  eq.  (2.2),  H  by  eq.  (2.3),  d^  by  eq.  (3.2),  and  S  by 
eq.  (3.3). 

Indices  (I,J)  of  one  grid  may  be  converted  to  indices  C I '  , J ' )  of  any  other 
grid  by 


1' 


(3.16) 


J' 


J' 

P 


(.1 


(3.17) 


Here  M  and  M1  are  the  mesh  factors  for  the  grids  with  points  (l,J)  and 
(l',J'),  respectively.  These  two  equations  apply  to  any  two  polar 
stereographic  grids  projected  on  the  same  hemisphere. 


3.1.2  POLAR  STEREOGRAPHIC  REFERENCE  GRIDS 


Polar  stereographic  reference  grids  exist  for  the  Northern  and  Southern 
Hemispheres  at  resolutions  to  provide  weather  support  from  the  synoptic  scale 
down  to  the  mesoscale.  For  a  polar  stereographic  projection  centered  on  a 
given  hemisphere,  the  domains  of  all  the  reference  grids  are  the  same;  they 
include  the  entire  hemisphere  on  which  the  grid  is  centered  and,  at  the 
corners  of  the  grids,  extend  well  beyond  the  Equator  into  the  opposite 
hemisphere.  (Fig.  3.6  provides  an  example  using  the  Northern  Hemispheric 
Whole-mesh  Reference  Grid.) 
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3.1. 2.1  THE  NORTHERN  HEMISPHERIC  WHOLE-MESH  REFERENCE  GRID 

The  Northern  Hemispheric  Whole-mesh  Reference  Grid  (sometimes  called  the 
Northern  Hemispheric  Whole-Mesh  Super  Grid)  is  a  65x65  whole-mesh  grid  based 
on  a  Northern  Hemispheric  polar  stereographic  map  with  standard  latitude 
60°N  (Fig.  3.6).  For  a  map  scale  of  1:30,000,000,  the  grid  can  be 
constructed  by  superimposing  a  65x65  array  of  points  on  a  map  with  one-half 
inch  between  points.  For  the  Northern  Hemispheric  Whole-mesh  Reference  Grid, 
eqs.  (3.1 )— (3 . 15)  apply  with  the  parameter  values  given  in  Table  3.1. 

The  (I,J)  indexing  convention  for  the  Northern  Hemispheric  Whole-mesh 
Reference  Grid  is  shown  in  Fig.  3.6.  The  column  index  I  varies  from  1  on  the 
left  of  the  grid  to  65  on  the  right,  whereas  the  row  index  J  varies  from  1  at 
the  top  of  the  grid  to  65  at  the  bottom.  Column  33  (I  =  33)  coincides  with 
the  100°E  and  80°W  meridians;  row  33  (J  =  33)  coincides  with  the  170°W 
and  10°E  meridians.  The  North  Pole  is  represented  by  point  (33,33). 

In  Fig.  3.6,  the  (whole-mesh)  grid-point  Cartesian  coordinate  system  (X,Y) 
defined  for  the  Northern  Hemispheric  polar  stereographic  projection  is 
superimposed  on  the  Northern  Hemispheric  Whole-mesh  Reference  Grid  with  the 
origin  of  the  system  at  the  North  Pole.  The  positive  X  axis  is  oriented  in 
the  direction  of  increasing  I  and  the  positive  Y  axis  in  the  direction  of 
dec reasing  J.  The  (X,Y)  coordinate  system  serves  as  the  frame  of  reference 
for  orientation  of  gridded  wind  components:  u  and  v  components  are  positive 
in  the  directions  of  increasing  X  and  Y,  respectively. 


Table  3.1.  Parameters  for  the  Northern  Hemispheric  Whole-mesh  Reference  Grid. 


PARAMETER 

DESCRIPTION 

H  =  +1 
(|)Q  =  60°N 
M  =  1 

de(at  <|’  =  60°N)  =  381  km 
1  1  I  ^  65 
1  <  J  <  65 


Northern  Hemispheric  projection 

standard  latitude 

whole-mesh  grid 

grid  length  at  60°N 

grid  points  in  the  I  direction 

grid  points  in  the  J  direction 

polar  location  (North  Pole) 


(Ip.Jp)  =  (33,33) 


3. 1.2.2  THE  SOUTHERN  HEMISPHERIC  WHOLE-MESH  REFERENCE  GRID 


The  Southern  Hemispheric  Whole-mesh  Reference  Grid  (also  called  the 
Southern  Hemispheric  Whole-mesh  Super  Grid)  is  a  65x63  whole-mesh  grid  based 
on  a  Southern  Hemispheric  polar  stereographic  map  true  at  60°S  (Fig.  3.7). 

For  a  map  scale  of  1:30,000,000,  the  grid  can  be  constructed  by  superimposing 
a  65x65  array  of  points  on  a  map  with  one-half  inch  between  points.  For  the 
Southern  Hemispheric  Whole-mesh  Reference  Grid,  eqs.  (3.1)— (3.13)  apply  with 
the  parameter  values  given  in  Table  3.2. 

The  (I,J)  indexing  convention  is  shown  in  Fig.  3.7.  1  extends  from  1  on 

the  left  of  the  grid  to  65  on  the  right  and  J  extends  from  1  on  the  top  of  the 
grid  to  65  at  the  bottom.  Column  33  (I  =  33)  coincides  with  the  100°E  and 
80°W  meridians;  row  33  (J  =  33)  coincides  with  the  170°W  and  10°E 
meridians.  The  South  Pole  is  located  at  grid  point  (33,33). 

In  Fig.  3.7,  the  (whole-mesh)  grid-point  Cartesian  coordinate  system  (X,Y) 
for  the  Southern  Hemispheric  polar  stereographic  projection  is  superimposed  on 
the  Whole-mesh  Reference  Grid  for  that  hemisphere.  The  origin  is  placed  at 
the  South  Pole  with  the  positive  X  axis  oriented  in  the  direction  of 
increasing  I  and  the  positive  Y  axis  oriented  in  the  direction  of  increasing 
J.  As  with  the  Northern  Hemispheric  grid,  this  coordinate  system  serves  as 
the  frame  of  reference  for  orientation  of  gridded  wind  components.  The  u 
component  is  positive  in  the  direction  of  increasing  X  and,  in  contrast  to  the 
Northern  Hemispheric  convention,  the  v  component  is  positive  in  the  direction 
of  decreasing  Y. 


Table  3.2.  Parameters  for  the  Southern  Hemispheric  Whole-mesh  Reference  Grid. 


PARAMETER 


DESCRIPTION 


H  =  -1 


Southern  Hemispheric  projection 


<f>0  =  60°S 


standard  latitude 


M  =  1 


dp( at  <|>  =  60°S)  =  381  km 


1  <  I  <  65 


whole-mesh  grid 

grid  length  at  60°S 

grid  points  in  the  I  direction 


1  <  J  <  65 


(ln,JD)  -  (33,33) 


grid  points  in  the  J  direction 
polar  location  (South  Pole) 
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J.i.2.3  FI  NKR-MESII  POLAR  ST  ERE  00  RAP  1 1 LC  REFERENCE  GRIDS 

Ln  representing  the  atmosphere  with  a  whole-mesh  grid,  approximation 
resuLts  because  sma l le r-sca Le  motions  cannot  be  correctly  represented.  In 
fact,  the  shortest  wavelength  resolvable  on  any  grid  is  equal  to  twice  the 
distance  between  its  grid  points.  Also,  wavelengths  shorter  than  about  four 
sriu  lengths  are  handled  poorly  by  numerical  prediction  models  (Haltiner  and 
Williams,  1980).  As  mentioned  previously,  the  whole-mesh  grid  length  is 
approximately  381  kin  (exactly  381  km  at  60°N  (or  °S))  on  the  earth's 
surface.  Consequently,  wavelengths  shorter  than  about  1600  km  are  poorly 
forecast  using  a  whole-me3h  grid  and  wavelengths  shorter  than  about  800  km  are 
not  represented  at  all. 

because  many  meteorological  applications  at  AFGWC  are  concerned  with 
scales  of  motion  much  smaller  than  either  1600  or  800  km,  f iner-resolution 
grids  iiad  to  be  developed.  With  the  whole-inesh  reference  grid  as  a  standard, 
a  family  of  finer-mesh  reference  grids  that  cover  the  same  total  area  as  the 
Whole-mesh  Reference  Grid  has  been  developed  for  both  the  Northern  and 
Southern  Hemispheres.  Figs.  3.8  and  3.9  for  the  Northern  and  Southern 
Hemispheres,  respectively,  include  reference  grids  at  half  mesh,  quarter  mesh, 
eighth  mesh,  and  sixty-fourth  mesh.  As  the  names  imply,  distance  between  grid 
points  on  these  finer-mesh  grids  is  defined  relative  to  the  reference 
whole-mesh  grid  length.  For  instance,  at  a  given  latitude,  distance  between 
adjacent  half-mesh  grid  points  is  one-half  the  distance  between  adjacent 
whole-mesh  points.  (For  example,  at  60°N  (or  °S),  the  whole-mesh  grid 
length  is  381  km,  so  that  the  half-mesh  grid  length  is  190.5  km).  The  AFGWC 
reference  grids  are  listed  by  mesh  size  (that  is,  whole-mesh,  half-mesh,  etc.) 
in  Table  3.3.  Also  included  in  the  table  are  the  mesh  factor  M,  which  is  the 
ratio  of  the  grid  length  to  the  whole-mesh  grid  length;  the  array  size,  where 
IxJ  is  the  total  number  of  grid  points;  (Ip,Jp),  which  are  the  (I,J) 
indices  of  the  pole;  and  the  distance  de  between  grid  points  given  in  both 
kilometers  and  nautical  miles.  The  basic  parameters  for  che  Finer-mesh  Polar 
Stereographic  Reference  Grids,  as  well  as  the  Whole-mesh  Polar  Stereographic 
Reference  Grids,  are  summarized  in  Table  3.4. 


Table  3.3. 

The  Polar  Stereogr3phic 

Reference  Grids. 

MESH  SIZE 

MESH  ARRAY  SIZE 

FACTOR  IxJ  =  TOTAL 

(M)  NUMBER 

OF  POINTS 

GRID  POINT 
LOCATION  OF 

THE  POLE 
dp>  Jp) 

DISTANCE  de  BETWEEN  GRID 
POINTS  AT  60°N  (OR  °S) 

KILOMETERS  NAUTICAL 

MILES* 

Who  le 

1 

65x65  =  4,225 

(33,33) 

381. 

205.72354 

Half 

1/2 

129x129  =  16,641 

(65,65) 

190.5 

102.86177 

Quart  e  r 

1/4 

257x257  =  66,049 

(129,129) 

95.25 

51.43088 

Eighth 

1/8 

513x513  =  263,169 

(257,257) 

47.625 

25. 71544 

Sixty-  fourtli 

1/64 

4097x4097  =  16,785,409 

(2049,2049) 

5.953125 

3.21443 

*based  on  the  international  nautical  mile  for  which  1  nm  =  1852  m. 
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Figure  3.8. 


,00,  — j  —  - 1— 1 — j— 

i  ..  ...  n  44-3^ 


\ww4#^  / 

\  :■■:  -  <>  /y  /VT^,-  \  \  \  v\  <X  x  x  :• 

\\  .  ‘  /  T  l  *  •  >4  ^Si  ’  \  „  \  '  «  .  '■Z 


/,  '  ‘  V,  /  .  .  \  \  r  '  \ 


\  •  V-  •  *  V-  J  /  •  ^-vjrii  -A  ;  '•  \ 

V/,  ,  •  •  A  .^r  ,  . « ^  • 

XV.  ■  /  Xf>.  '  V/>  3- 


/  ;  i: 


The  (I,J)  indices  for  the  Whole-mesh  and  Finer-mesh 
Reference  Grids  of  the  Northern  Hemisphere.  The  subscripts 
denote  whole-mesh  (1W,JW),  half-mesh  (l2»J2^> 
quarter-mesh  (l^J^),  eighth-mesh  (Lg.Jy),  and 
sixty-fourth-mesh  (l6P,Jfc>q)  indices.  The  who le-nwdi 
grid  spacing  is  displayed  in  the  upper-left  corner. 
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figure  3.9. 


The  (I,J)  indices  for  the  Whole-mesh  and  Finer-mesh 
Reference  Grids  of  the  Southern  Hemisphere.  The  subscripts 
denote  whole-mesh  (IW,JW),  half-mesh  (l2,J2)> 
quarter-mesh  (14^4)  >  eighth-mesh  (Is.Js) ,  and 
sixty- fourth-mesh  (I64,J64)  indices.  The  whole-mesh 
grid  spacing  is  displayed  in  the  upper-left  corner. 


Table  3.4.  Parameters  for  all  the  Polar  Stereographic  Reference  Grids. 


PARAMETER 

DESCRIPTION 

H 

hemisphere  indicator  (see  eq.  2.3) 

()>  =  (H)(60°  latitude) 

0 

standard  lutitude 

M 

mesh  factor  (see  Table  3.3) 

de  (at  <f>  =  (H)(60°  latitude)) 

distance  between  grid  points  on  the 
earth's  surface  at  (11)  (60°  latitude) 

( see  eq.  3.4) 

I 

grid-point  index  in  the  I  direction 
(see  Table  3.3) 

J 

grid-point  index  in  the  J  direction 
( see  Table  3.3) 

(I p  >  0  p ) 

grid-point  location  of  the  pole 
(see  Table  3.3) 

The  (I,J)  indexing  convention  and  the  geographical  orientations  of  the  two 
reference  grids  also  apply  to  their  finer-mesh  counterparts.  The  I  index  increases 
from  left  to  right  across  the  grid  and  the  J  index  increases  from  the  top  of  the 
grid  to  the  bottom.  If  dnp,Jnp)  represents  the  (I,J)  indices  of  the  North 
Pole,  column  1  =  I  coincides  with  the  100°E  meridian  in  the  top  half  of  the 
grid  and  with  the  80°W  meridian  in  the  bottom  half  of  the  grid;  row 
coincides  with  the  170°W  meridian  in  the  left  half  and  with  the  10°E  meridian 
in  the  right  half.  If  (I8p,Jsp)  represents  the  (I,J)  indices  of  the  South 
Pole,  column  I  =  Igp  coincides  with  80°W  in  the  top  half  of  the  grid  and  with 
100°E  in  the  bottom  half  of  the  grid;  row  J  =  Jgp  coincides  with  170°W  on  the 
left  and  with  10°E  on  the  right.  Indices  (I,J)  of  one  grid  may  be  converted  to 
indices  of  another  grid  by  eqs.  (3.16)  and  (3.17). 

The  orientations  of  the  grid-point  Cartesian  coordinate  systems  (X,Y)  for 
the  Finer-mesh  Reference  Grids  of  the  Northern  and  Southern  Hemispheres  are 
the  same  as  for  the  Whole-mesh  Reference  Grids  (see  Figs.  3.6  and  3.7). 

Equations  for  (X,Y)  coordinates  in  multiples  of  the  appropriate  grid  length 
d^  are  obtained  from  eqs.  (3.7)  and  (3.8)  with  H  =  +1  and  -  1  for  the 
Northern  and  Southern  Hemispheres,  respectively.  Gridded  wind  components  are 
oriented  using  the  same  conventions  as  the  Whole-mesh  Reference  Grids. 

Latitude  and  longitude  may  be  determined  from  a  grid  point's  (1,J) 
coordinates  from  eqs.  (3.11)  and  (3.12)  with  the  mesh  factor  M  specified  in 
Table  3.3.  The  (I,J)  coordinates  of  any  point  on  the  earth's  surface  can  be 
determined  by  eqs.  ( 3 . 1 3 )- ( 3 . 15 ) . 
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In  summary,  eqs.  (3. 1)— (3.  15)  apply  to  the  Finer-mesh  Reference  Grids,  as 
well  as  the  other  polar  stereographic  grids.  Eqs.  (3.16)  and  (3.17)  provide  a 
means  of  relating  indices  of  different  grids. 


1.1.  i  POLAR  STEREOGRAPHIC  GRIDS  FOR  NUMERICAL  WEATHER  ANALYSIS  AND  FORECASTING 

A  number  of  grids  for  specific  applications  have  been  defined  for  various 
numerical  weather  analysis  and  forecasting  models  used  at  AFGWC.  These  grids 
are  areal  subsets  of  the  reference  grids  and  vary  from  small  window  (about 
3000  km  by  3000  km)  to  hemispheric  in  areal  coverage.  A  brief  overview  of 
each  applications  grid  will  be  presented  in  this  section.  Descriptions  of  the 
functions  of  the  models  are  presented  by  Tarbell  and  Hoke  (1979). 

Eqs.  (J.l)-(3.17)  apply  to  these  polar  st  e  reographi.c  grids  with  the  basic 
grid  parameters  summarized  in  Table  3.5. 

3.1. 3.1.  NORTHERN  AND  SOUTHERN  HEMISPHERIC  OCTAGONS 

both  whole-mesh  and  half-mesh  versions  of  Northern  and  Southern 
Hemispheric  octagonal  grids  are  used  operationally  at  AFGWC.  These  grids  were 
developed  to  meet  the  demands  for  higher-accuracy  forecasts  while  satisfying 
computer  constraints. 

3. 1.3. 1.1.  WHOLE-MESH  OCTAGONS 

The  Whole-mesh  Octagon  was  the  first  grid  used  for  numerical  weather 
prediction  at  AFGWC.  Octagon  grids  have  been  defined  for  both  the  Northern 
and  Southern  Hemispheres  (Figs.  3.10  and  3.11).  The  grid  size  is  97x51  with 
the  corner  portions  of  the  rectangle  deleted.  That  is,  for  the  first  row  (J  = 
1)  of  the  Octagon,  only  the  grid  points  (15,1)  through  (33,1)  are  included;  in 
row  two,  points  (19,2)  through  (39,2)  are  used.  Similarly,  each  successive 
row  includes  two  more  grid  points  until  a  full  complement  of  97  grid  points  is 
attained  in  row  15.  After  row  37,  in  which  points  (1,37)  through  (97,37)  are 
included,  each  successive  row  deletes  two  points  (one  at  each  end).  Row  38 
includes  (2,38)  through  (96,38),  row  39  includes  (3,39)  through  (95,39),  and 
so  on.  The  final  row  of  the  Octagon  includes  points  (15,51)  through  (33,51). 

Information  is  actually  stored  for  the  octagon  on  a  rectangular  grid.  The 
information  in  the  corners  of  this  grid,  outside  the  octagon,  is  often 
meaningless  and  may  simply  be  assigned  a  constant  value  or  may  be  extraneous 
values  left  over  from  prior  computations.  For  example,  for  the  Northern 
Hemispheric  whole-mesh  analysis  data  base  that  is  designated  at  AFGWC  as 
literal  RTGWCA,  data  in  the  corners  of  the  grid  should  usually  be  ignored. 

The  earth's  pole  is  represented  by  grid  point  (29,26)  for  both  the 
Northern  and  Southern  Hemispheric  Whole-mesh  Octagons,  whereas  on  the 
Whole-mesh  Reference  Grid  the  pole  location  is  (33,33).  Whole-mesh  Octagon 
indices  may  be  converted  to  their  Whole-mesh  Reference-grid  values  by  eqs. 
(3.16)  and  (3.17),  where  (I,J)  are  the  indices  on  the  Whole-mesh  Octagon, 
(I'jJ1)  the  indices  on  the  Whole-mesh  Reference  Grid,  and  M  =  M1  =  1  being  the 
respective  mesh  factors. 

Various  parameters  of  the  Whole-mesh  Octagons  are  summarized  in  Table  3.5. 
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3.  1.3. 1.2.  HALF-MESH  OCTAGONS 

The  Half-mesh  Octagons  (Figs.  3.10  and  3.11),  which  contain  93x10)  points, 
are  subsets  of  the  Half-mesh  Reference  Grids.  As  with  the  Whole-inesh 
Octagons,  grid  points  in  the  corner  portions  of  the  rectangle  are  not  used. 
Half-mesh  Octagon  indices  (I,J)  may  be  converted  to  their  Half-mesh 
Reference-grid  values  (I  *  ,  J ' )  using  eqs.  (3.16)  and  (3.17),  where  Op,Jp) 

=  (47,51),  (Ip'jjp1)  =  (65,65),  and  M  =  M'  =  1/2.  Whole-mc3h  Octagon 
indices  (I,J)  may  be  converted  to  their  Half-mesh  Octagon  values  (I * , J 1 )  also 
using  eqs.  (3.16)  and  (3.17)  with  (Ip,Jp)  =  (24,26),  (Ip'.Jp1)  = 

(47,51),  M  =  l,  and  M*  =  1/2. 

Various  parameters  of  the  Half-mesh  Octagons  are  listed  in  Table  3.5. 


3. 1.3. 1.3  HISTORICAL  NOTE 

Size,  shape,  and  other  attributes  of  the  Octagons  were  not  arbitrarily 
specified.  Development  of  the  grid  system  began  in  the  mid-1950s  at  the  Joint 
Numerical  Weather  Prediction  Unit  ( JNWPU) ,  Suitland,  Maryland  (H.  A.  Bedient, 
personal  communication).  Many  of  these  attributes  resulted  from  constraints 
imposed  by  computer  systems  in  those  early  days  of  operational  numerical 
weather  prediction.  The  line  printer  played  an  important  role,  for  example. 
Because  this  device  printed  characters  at  a  density  of  10  characters  per  inch 
horizontally  and  vertically  on  the  page,  computer-generated  gridded  output  was 
conveniently  plotted  using  a  grid  spacing  of  1/2  inch  for  the  Whole-mesh 
Octagon.  On  a  polar  stereographic  projection  with  map  scale  M  of  1:30,000,000 


Table  3.5.  Hemispheric  Grids  for  Numerical  Weather  Analysis  and  Forecasting. 
The  pole  location  is  with  respect  to  each  particular  grid,  not  the  Reference 
Grids.  NH  indicates  Northern  Hemisphere  (H  =  +1,  >)'()  =  60°N)  and  SH 
indicates  Southern  Hemisphere  (H  -  -1,  <l»  -  60°S). 


GRID 

DOMAIN 

MESH 

FACTOR 

(M) 

ARRAY 

IxJ  = 

NUMBER 

SIZE 

;  TOTAL 

OF  POINTS 

GRID  POINT 

LOCATION  OF 

THE  POLE 

(lp,Jp) 

Whole-mesh  Octagon 

NH/SH 

1 

47x51  = 

2,397 

(  24 , 26  ) 

Half-mesh  Octagon 

NH/  SH 

1/2 

93x101  = 

9,393 

(47, 51) 

TRONEW 

NH/SH 

1/2 

128x128  = 

16,384 

(65,65) 

AWSPE 

NH  only 

1 

53x57  = 

3,021 

(27,29) 

SIXLVL 

NH/SH 

1 

51x51  = 

2,601 

( 26 , 26  ) 

3DNEPH 

NH/SH 

1/8 

512x512  = 

262,144 

(257,257) 

SGDB 

NH/SH 

1/64 

4096x4096  = 

=  16,777,216 

(2049,2049) 
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F igure  3.11. 


The  Whole-mesh  and  Half-mesh  Octagons  for  Lhe  Southern 
Hemisphere.  The  domain  of  the  Reference  Grids  (solid  outer 
border)  is  also  plotted.  The  indices  Ow,Jw)  ami 
( 1 2  ,  J 2  )  designate  the  coordinates  for  the  whole-mesh 
and  half-mesh  versions,  respectively,  with  respect  to  the 
Octagons  and  the  Reference  Grids.  The  dashed  lines 
indicate  that  data  on  the  Octagons  are  actually  stored  in  a 
rectangular  data  base.  The  whole-mesh  grid  spacing  is 
displayed  in  Lhe  upper-left  corner. 
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and  true  at  60°N ,  which  was  a  common  projection  then  and  now,  1/2  inch 
corresponds  to  exactly  381  kin  on  the  earth  at  60°N . 

Because  of  the  severe  computer  memory  constraints  of  the  IBM  704  computer 
system  used  at  that  time,  the  number  of  grid  points  had  to  be  less  than  2048. 
The  1477  grid  points  of  the  Whole-mesh  Octagon  fit  this  requirement.  The 
octagon  was  selected  instead  of  a  more  circular  polygon  to  simplify 
development  of  the  numerical  models.  The  47x51  octagon  was  selected  instead 
of  an  equilateral  octagon  in  part,  once  again,  because  of  the  line  printer; 
the  47  grid  points  in  the  I  direction  (left  to  right)  fit  on  two  computer 
pages  side  by  side.  (At  that  time,  only  120  characters  could  be  printed  per 
line  on  a  page.)  Although  up  to  48  grid  points  would  fit  on  two  pages,  this 
larger  size  was  not  used  because  of  the  desire  to  have  the  North  Pole  located 
at  a  grid  point,  thereby  necessitating  a  grid  with  odd  dimensions.  The 
maximum  odd  number  of  grid  points  that  can  be  included  in  the  J  direction  (top 
to  bottom)  without  exceeding  the  2048  limit  for  total  points  is  51. 

The  standard  longitude  A0  was  selected  as  10°E  to  maximize  the  number  of 
important  data  points  near  the  boundaries  that  could  be  included  within  the 
domain  of  the  Northern  Hemispheric  Whole-mesh  Octagon.  For  example,  a 
rotation  of  the  Octagon  by  less  than  15°  longitude  to  the  west  places  Dakar 
outside  the  domain,  while  a  rotation  of  about  25°  longitude  to  the  west 
eliminates  Manila. 

The  Southern  Hemispheric  Whole-mesh  Octagon  was  then  set  up  to  be 
analogous  to  its  Northern  Hemispheric  counterpart,  and  the  Half-mesh  Octagons 
were  devised  to  meet  the  demands  for  higher  resolution. 

3.  1.3.2.  AWSPE  AND  SIXLVL  HEMISPHERIC  GRIDS 

The  grid  for  storing  information  pertinent  to  the  Air  Weather  Service 
Primitive  Equation  Model  (AWSPE)  is  a  53x57  subset  of  the  Whoie-mesh  Reference 
Grid.  This  subset  grid  is  currently  defined  and  used  operationally  for  the 
Northern  Hemisphere  only  (Fig.  3.12).  Although  the  (I,J)  indexing  convention 
of  the  grid  for  storing  AWSPE  output  is  as  shown  in  Fig.  3.12,  in  the  actual 
indexing  convention  interna  1  to  the  AWSPE  computer  program  the  lower-left 
corner  is  grid  point  (1,1). 

The  grid  for  the  Six-level  Baroclinic  Prediction  Model  (SIXLVL)  is  a  51x51 
subset  of  the  Whole-mesh  Reference  Grid.  SIXLVL  grids  for  both  the  Northern 
and  Southern  Hemisphere  (Figs.  3.13  and  3.14)  are  used  operationally. 

All  conventions  and  equations  described  in  the  section  on  reference  grids 
apply  for  the  AWSPE  and  SIXLVL  grids.  Indices  of  grid  points  for  these  grids 
may  be  converted  to  their  Whole-mesh  Reference-grid  values  using  information 
from  Tables  3.3  and  3.5  and  eqs.  (3.16)  and  (3.17). 

3.  1.3.3.  3DNEPH  GRIDS 

The  grids  for  the  Three-dimensional  Nephanalysis  Model  (3DNEPH)  are 
512x512  subsets  of  the  513x513  Eighth-mesh  Reference  Grids.  The  bottom  row 
(J  =  513)  and  the  right-most  column  (I  =  513)  are  not  included  in  the 
hemispheric  subset  grids  used  by  the  3DNEPH  (Figs.  3.15  and  3.16). 
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Figure  3.12. 


The  grid  (inner  solid  lines)  for  the  Air  Weather  Service 
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Figure  3.13. 


The  grid  (inner  .solid  lines)  lor  the  Six-lev 
Prediction  Model  (SIXLV1.)  for  L he  Northern  H 
domains  of  the  Northern  Hemispheric  Referenc 
outer  border),  the  AWSi’K  grid  (outer  dashed 
Octagon  grids  (inner  dashed  lines)  are  also 
indices  (IW,JW)  designate  the  coordinates  fo 
Whole-mesh  Reference  Grid  and  the  SIXLVI.  gn 
whole-meali  grid  spacing  ia  displayed  in  the 
cor  lie  r . 
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F  igure  3.15. 
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The  grid  system  ot  tin-  <  ree-d  uiiea;!  iona  1  Nephanalysis 
(3DNEPH)  for  the  Nortiu.  in  ilemi  spin:  re .  The  (l,J)  indexing 
convention  is  the  same  as  that  of  the  Eighth-mesh  Reference 
Grid.  The  indices  (lw<  5W1  and  (Ig.Jg)  denote  the 
whole-mesh  and  eighth-mesh  coordinates,  respectively,  for 
the  Reference  Grids  and  the  3DNEPH.  The  64  3DNEPH  boxes 
are  indicated;  the  four  corner  boxes  are  not  used.  The 
(ln,jn)  coordinate  axes  in  tin  lower-right  corner  show 
the  3DNEPH  Box  Ij  convent  ion  for  .i  sample  3DNEPH  box.  Each 
3DNEPH  box  contains  the  grid  points  on  the  upper  and  left 
boundaries  of  the  box,  but  not  th  lower  and  right 
boundaries.  The  labels  (.64.H/5)  and  (512)  on  the  axes 
respectively  indicate  the  whole-mesh  and  eighth-mesh  values 
of  the  last  row  of  the  3DNFPH  grid  system,  which  has  one 
less  row  than  the  Eighth-mesh  Reference  Grid.  The 
whole-mesh  grid  spacing  is  displayed  in  the  upper-left 
corne  r . 


Figure 


The  grid  system  of  the  Three-dimensional  Nephanalysis 
(3DNEPH)  for  the  Southern  Hemisphere.  The  (l,J)  indexing 
convention  is  the  same  as  that  of  the  Eighth-mesh  Reference 
Grid.  The  indices  (IW,JW)  and  (Ig-Jg)  ^note  the 
whole-mesh  and  eighth-mesh  coordinates,  respectively,  tor 
the  Reference  Grids  and  the  3DNEP11.  The  64  JDNEPH  boxes 
are  indicated;  the  four  corner  boxes  are  not  used.  The 
( T.  j_)  coordinate  axes  in  the  lower-right  corner  show 
the  3DNEPH  Box  IJ  convention  for  a  sample  3DNEPH  box.  Each 
3DNEPH  box  contains  the  grid  points  on  the  upper  and  It  t 
boundaries  of  the  box,  hut  not  the  lower  and  right 
boundaries.  The  labels  (64.873)  and  (312)  on  the  axes 
respectively  indicate  the  whole-mesh  and  eighth-mesh  values 
of  the  last  row  of  the  3DNEPH  grid  system,  which  has  one 
less  row  than  the  Eighth-mesh  Reference  Grid.  The 
whole-mesh  grid  spacing  is  displayed  in  the  upper  e  l 


corner. 


Information  stored  on  3DNEPH  grids  is  assumed  to  be  representative  of 
conditions  surrounding  each  grid  point.  This  is  in  contrast  to  the  convention 
for  the  Satellite  Global  Data  Base  (see  Section  3. 1.3.5). 


As  noted  by  Fye  (1978),  the  array  size  (262,14a  points)  of  the  3DNEPH 
hemispheric  grid  "requires  partitioning  to  facilitate  automated  processing." 
For  each  hemisphere,  the  512x512  array  is  therefore  divided  into  64  areas 
called  3DNEPH  boxes  (Figs.  3.15  and  3.16),  denoted  here  by  Bn.  The  four 
corner  boxes,  which  lie  off  the  hemisphere,  are  not  used.  The  box  number  Bn 
can  be  calculated  from  indices  (Iy,Jg)  of  Lhe  Eighth-mesh  Reference  Grid 
according  to 

b  =  1  +  1  NT ( 1,-1,  6)0  +6  [  TNT( JQ-1  ,  6'Oi  (3.18) 

1 1  Cj  o 

tor  1  Tig-  512  and  1  £  Jg  1_  512.  INT  designates  a  function  that 
provides  the  integer  portion  of  the  division  of  the  first  argument  by  the 
second.  Also,  we  let  MOD  designate  a  mathematical  function  in  which  the 
functional  value  is  the  remainder  of  the  division  of  the  first  argument  by  the 
second  argument.  For  example,  MOD(7,2)  =  1  and  INT(7,2)  =  3;  MOD(57,8)  =  1 
and  INT(57,8)  =  7.  Note  that  1  .*1  Bn  0  64. 

Within  each  box  an  indexing  convention  (In,Jn),  called  the  3DNEPH  Box 
1J  convention,  is  used  to  designate  the  64x64  array  of  grid  points  (see  Fig. 
3.17).  Box  IJ  indices  (In,Jn)  are  related  to  their  corresponding 
Eighth-mesh  Reference-grid  indices  (Ig,Jg)  by 

lfl  =  I'^-l  ,  8)  i  ,  (3.19) 


i8  -  .Jn  +  6M  lNT(Hn-l,  8)]  ,  (3.20) 

where  Bn  is  the  3DNEPH  box  number,  1  ^In  <  64,  and  1  ^  Jn  0&4. 

Eighth-mesh  Reference  Grid  indices  may  be  transformed  to  3DNEPH  Box  IJ  indices 
for  box  Bn  by 

1  =  1  +  MOD ( 1 - 1  ,  bh)  ,  (3.21) 

n  o 


Jit  =  1  +  MOD  (<J - 1  ,  64) 


(3.22) 


3.  1.3.4.  TRONEW  HEMISPHERIC  GRIDS 

The  grids  for  the  Tropical  Cloud  Forecast  Model  (TRONEW)  are  the  same  as 
the  129x129  Half-mesh  Reference  Grids  (Figs.  3.18  and  3.19).  Although  the 
TRONEW  grids  are  129x129,  useful  information  is  only  stored  on  the  128x128 
subsets.  The  bottom  row  (J  =  129)  and  the  right-most  column  (I  =  129)  of  the 
llalf-mesh  Reference  Grids  are  not  included  in  the  Northern  and  Southern 
Hemispheric  grids  used  by  TRONEW.  All  conventions  and  equations  described  in 
the  sections  on  the  generalized  equations  (Section  3.1.1)  and  the  finer-mesh 
polar  slereographic  reference  grids  (Section  3. 1.2.3)  apply  for  the  TRONEW 
grids.  Various  parameters  of  these  grids  are  listed  in  Table  3.5. 
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Figure  3.17. 
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The  relationship  of  3DNEPH  Box  IJ  indices  (ln>Jn)  to 
Whole-mesh  Reference  Grid  indices  Ow,Jw)  and 
Eighth-mesh  Reference  Grid  indices  (lg,Jg)  for  an 
example  3DNEPH  box  (box  number  Bn  =  4 4).  Note  that  grid 
points  on  the  lower  and  right  boundaries  (solid  lines)  of 
the  box  are  located  in  adjacent  boxes.  The  dashed  lines 
indicate  the  actual  lower  and  right  boundary  rows  of  the 
box,  with  (I,J)  indices  given  in  parentheses  on  t  lie  axes. 


iwMsaBaea^:-  ~ 


The  Northern  Hemispheric  grid  of  the  Tropical  Cloud 
Forecast  Model  (TRONEW).  The  grid  is  half  mesh.  The  (I,J) 
indexing  convention,  designated  ( 1 2  ,  J  2  ^ »  i-s  the  same  as 
that  of  the  Half-mesh  Reference  Grid  (solid  outer  border). 
The  dashed  lines,  labeled  (128)  on  the  axes,  indicate  that 
meaningful  data  are  stored  on  a  128x128  domain,  although 
the  size  of  the  data  base  is  129x129.  The  whole-mesh  grid 
spacing  is  displayed  in  the  upper-left  corner. 
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Figure 


3.19.  The  Southern  Hemispheric  grid  of  the  Tropical  Cloud 

Forecast  Model  (TRONEW) .  The  grid  is  half  mesh.  The  (l,j) 
indexing  convention,  designated  ( 1 2  , J 2 ^  >  is  t'ie  same  as 
that  of  the  Half-mesh  Reference  Grid  (solid  outer  border). 
The  dashed  lines,  labeled  (128)  on  the  axes,  indicate  that 
meaningful  data  are  stored  on  a  128x128  domain,  although 
the  size  of  the  data  base  is  129x129.  The  whole- mesh  grid 
spacing  is  displayed  in  the  upper-left  corner. 
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1.  1.3.5  SATELLITE  GLOBAL  DATA  BASE  GRIDS 


The  Satellite  Global  Data  Base  (SGDB)  uses  a  4096x4096  subset  of  the 
Sixty- fourth-mesh  Reference  Grid.  The  bottom  row  (J  =  4097)  and  the 
right-most  column  (I  =  4097)  of  this  reference  grid  are  not  included  in  the 
two  hemispheric  subsets  of  the  SGDB  (Figs.  3.20  and  3.21). 

Like  the  3DNEPH  grids,  the  SGDB  hemispheric  grids  are  partitioned  into 
boxes  to  facilitate  computer  processing  of  satellite  data.  For  each 
hemisphere,  there  are  4096  SGDB  boxes  (Figs.  3.20  and  3.21),  each  consisting 
of  a  64x64  array  of  points,  for  a  total  of  16,777,216  grid  points  per 
hemisphere.  Each  SGDB  box  contains  one  Whole-mesh  Reference  Grid  Point 
(IW,JW),  which  coincides  with  the  s ixty- fourth-mesh  point  in  the 
upper-left  corner  of  the  SGDB  box  (Fig.  3.22).  Imagery  is  not  stored  in  those 
boxes  that  lie  off  the  hemisphere.  SGDB  boxes  are  numbered  according  to 


I!  =  (1  -1)  +  64  (.!  -1) 

H  W  W 


11  64 

—  w 

1  .J  •  64 
w 


(3.23) 


where  (IW,JW)  are  the  indices  of  the  Whole-mesh  Reference  Grid  Point  at 
the  SGDB  point  in  the  upper-left  corner  of  the  box  and  Bg  is  the  SGDB  box 
number.  The  box  number  can  be  alternately  calculated  from  the 
Sixty-fourth-mesh  Reference  Grid  Indices  (I64.J64)  according  to 

B„  =  TNT(  J  6U- 1 ,  64  )  +  6b  [TNT(.lf)j-l  ,  64)]  (3.24) 

for  1  <_  1(j4  <_  4096  and  1  <_  J^4  4096.  Note  that  0  ^  Bs  5.  4093. 


Within  each  box  an  indexing  convention,  called  the  SGDB  Box  IJ  convention, 
is  used.  SGDB  Box  IJ  indices  (Ig,Js)  are  related  to  their 
Sixty- fourth-mesh  Reference  Grid  indices  (1-64, J64)  by 


1 64  =  l,.  +  64 1  MOD  (  ,  64  )j 


(3.25) 


J64  =  J,.  +  64|INT(B_,  64)] 


(3.26) 


where  1  <_  Ig  ^64,  1  Js  64,  and  MOD  and  INT  are  defined  in  Section 
1.1. 3. 3.  Si xty- fourth-mesh  Reference  Grid  indices  may  be  transformed  to  SGDB 
Box  IJ  indices  for  box  Bg  by 


+  Mi  Jilt.  I  -  1  ,  04  ) 

(  Wj 


i  +  MonU^-i,  64) 


(3.27) 

(3.28) 
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The  grid  system  for  the  Satellite  Global  Data  Base  (SGDB) 
for  the  Northern  Hemisphere.  The  (l,J)  indexing  convention 
is  the  same  as  that  of  the  Sixty-fourth-mesh  Reference 
Grid.  The  indices  (lw,Jv)  and  (l64> J64)  denote  the 
coordinates  of  the  Whole-mesh  Reference  Grid  and  the 
Sixty-fourth-mesh  Reference  Grid,  respectively,  whereas  the 
indices  ( Iw* ,JW’ )  and  (164 ' *  Jfe4  ' )  denote  the 
numbering  convention  used  in  the  satellite  data  processing 
software  (SYNAPS).  Examples  of  the  4096  SGDB  boxes,  which 
are  one  whole-mesh  grid  length  on  a  side,  are  given  in  the 
upper-left  corner  of  the  domain.  Imagery  is  not  stored  in 
those  boxes  that  lie  off  the  hemisphere.  Each  SGDB  box 
contains  the  rows  of  grid  points  on  the  upper  and  Left 
boundaries  of  the  box,  but  not  the  lower  and  right 
boundaries.  The  index  values  in  parentheses  on  the  axes 
indicate  the  values  of  the  last  row  of  the  SGDB  grid 
system,  which  has  one  less  row  than  the  Sixty- fourtli-mesli 
Reference  Grid.  Over  the  U.S.,  box  number  2977  (Bs  = 

2977)  is  given  as  an  example  and  is  described  in  more 
detail  in  Fig.  3.22 . 
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Figure  3.21. 


\y  ^  -  *  . 

■S'1..-  - 

>i-K.v  ,:^:0 

i'j,  '  '■  vl  *  * 

■  sp-V'0sa*,<"' 

Jv-*-*" 


,  ••  • 


Jp  X  v 


v  v  v  sfi  ■  ■ 

v  •  t'vii 

\  £.;a  ■  ai 

V«>r/  A:  • ..,:  J 


,  a 

\ _  '■•  i.'  *■■'»' 

r-V  , 

.  •  t  '" 

>  .  '■> 

>  ,•  /'••>•. . 


'YV  , 


■r*-**' 


>.u\  .  '  ■‘•-  v 

x  •*>  *.  «. 


The  grid  system  for  the  Satellite  Global  Data  Base  (SGDB) 
for  the  Southern  Hemisphere.  The  (I,J)  indexing  convention 
is  the  same  as  that  of  the  Sixty- fourth-mesh  Reference 
Grid.  The  indices  dw,Jw)  and  (l64> J64)  denote  the 
coordinates  of  the  Whole-mesh  Reference  Grid  and  the 
S lxty- fourth-mesh  Reference  Grid,  respectively,  whereas  the 
indices  ( 1 w 1 , Jw '  )  and  C 1 6 4 1 , 3&4  ’  ^  denote  the 
numbering  convention  used  in  the  satellite  data  processing 
software  (SYNAPS).  Examples  of  the  4096  SGDB  boxes,  which 
are  one  whole-mesh  grid  length  on  a  side,  are  given  in  the 
upper-left  corner  of  the  domain.  Imagery  is  not  stored  in 
those  boxes  that  lie  off  the  hemisphere.  Each  SGDB  box 
contains  the  rows  of  grid  points  on  the  upper  and  left 
boundaries  of  the  box,  but  not  the  lower  and  right 
boundaries.  The  index  values  in  parentheses  on  the  axes 
indicate  the  values  of  the  last  row  of  the  SGDB  grid 
system,  which  has  one  less  row  than  the  Sixty-fourth-mesh 
Reference  Grid. 


The  relationship  of  SGDB  Box  IJ  indices  (1Q,J8)  to 
Whole-mesh  Reference  Grid  indices  (IW,JW)  and 
Sixty-fourth-mesh  Reference  Grid  indices  ( If,/, ,  Jf,/, )  f< 
an  example  SGDB  box  (box  number  Bg  =  2977).  The  ind 
ds'.Jg'),  (IW',JW'),  and  (l64,»J64')  denote 
the  numbering  convention  used  in  the  satellite  data 
processing  software  (SYNAPS).  Grid  points  on  the  lo- 
right  boundaries  (solid  lines)  of  the  box  are  locate- 
adjacent  boxes.  The  dashed  lines  indicate  the  actua 
and  right  boundary  rows  of  the  box,  with  (I,J)  indie 
given  in  parentheses  on  the  axes. 
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Satellite  data  processing  software  (SYNAPS)  designates  the  origin  of  the 
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1=1-1 
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where  (IW',JWI),  ( 1^4 '  ,  '  )  ,  and  (Ig'.Jg1)  are  whole-mesh 

indices,  sixty-fourth-mesh  indices,  and  SGDB  Box  IJ  indices,  respectively,  as 
used  in  SYNAPS. 

Information  stored  on  SGDB  grids  is  not  representative  of  conditions 
surrounding  each  grid  point,  as  i3  the  case  for  other  AFGWC  grids,  such  as  the 
3DNEPH  grids.  Instead,  SGDB  information  stored  at  a  Sixty-fourth-mesh 
Reference  Grid  Point  (1^4,354)  represents  the  conditions  in  the  area 
bounded  by  that  grid  point  and  ( I(,4+ 1 ,  J64) ,  (l64>J64+D,  and  (I64+1»J64+D* 

3.  1.3.6.  STATIONARY  WINDOW  GRIDS. 

Smaller-area  subsets,  or  window  grids,  are  also  used  in  AFGWC  models. 

These  window  grids  include  stationary  windows  for  the  North  American, 

European,  and  Asian  continents  and  movable  windows  of  nearly  constant 
areal  extent.  Each  window  is  defined  relative  to  a  specific  reference  grid 
and  obeys  all  conventions  described  for  that  reference  grid.  Table  3.6 
details  important  parameters  of  each  stationary  window  grid.  Locations  of  the 
corner  points  with  respect  to  the  coordinates  of  the  Whole-mesh  Octagon  were 
derived  using  eqs.  (3.16)  and  (3.17)  with  (I,J)  being  the  coordinates  of  the 
points  with  respect  to  each  particular  window  grid,  (I  * ,  J ' )  the  coordinates  on 
the  Whole-mesh  Octagon,  M  =  1/2  for  half-mesh  grids,  M'  =  1  for  the  Whole-mesh 
Octagon,  (Ip'.Jp')  =  (47,51)  from  Table  3.5,  and  (Ip,Jp)  from  Table 
3.6  for  the  respective  stationary  window  grids. 


3. 1.3.6. I  STATIONARY  WINDOW  GRIDS  FOR  THE  HALF-MESH  ANALYSIS  MODELS. 


The  three  windows  used  for  the  Half-mesh  Analysis  Models  are  shown  in  Fig. 
3.23  (solid  lines)  and  summarized  in  Table  3.6.  These  grids  are  defined 
relative  to  the  Half-mesh  Reference  Grid  and  include  a  37x39  half-mesh  window 
for  North  America  and  35x41  half-mesh  windows  for  the  European  and  Asian 
areas.  The  (I,J)  coordinate  axes  are  as  defined  for  the  Half-mesh  Reference 
Grid. 
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Table  3.6.  Stationary  Window  Grids.  The  grid- point  location  of  the  pole  is 
with  respect  to  each  particular  stationary  window  grid,  not  the  Reference 
Grids.  (The  pole  is  only  within  the  domain  of  the  European  and  Asian 
Half-mesh  Analysis  Model  grids.)  The  mesh  factor  M  for  all  these  grids  Is 
1/2.  The  (I,J)  coordinates  (on  the  whole-mesh  octagon)  of  the  corners  of  the 
window  grids  can  be  determined  using  other  information  in  this  table  and  eqs. 
(3.16)  and  (3.17). 


DOMAIN 

.ARRAY  SIZE 

IxJ  =  TOTAL 
NUMBER  OF 

POINTS 

GRID  POINT 

LOCATION  OF 

THE  POLE 
(  i  P  >  dp) 

(I,J)  COORDINATES  (OF  THE 
WHOLE- MESH  OCTAGON)  FOR  THE 
CORNERS  OF  THE  WINDOW  GRIDS 

Corner  Points 

Upper  Upper  Lower  Lower 
Left  Right  Left  Right 

Stationary  Window  Grids  for  Half-mesh  Analysis 

Models 

N.  America 

37x39  1443 

(27,-3) 

(11,28) 

(29,28) 

(11,47) 

(29,47) 

Europe 

35x41  =  1435 

(  1,29) 

(24,12) 

(41,12) 

(24,32) 

(41,32) 

Asia 

35x41  =  1435 

(25,41) 

(12,  6) 

(29,  6) 

(12,26) 

(29,26) 

Boundary-layer  Model  Grids 

U.S. 

29x27  =  783 

(23,-13) 

(13,33) 

(27,33) 

(13,46) 

(27,46) 

Europe 

29x35  =  1015 

(-3,27) 

(26,13) 

(40,13) 

(26, 30) 

(40,30) 

Asia 

29x35  =  1015 

(23,39) 

(13,  7) 

(27,  7) 

(13,24) 

(27,24) 

3. 1.3. 6. 2  BOUNDARY-LAYER  MODEL  GRIDS. 

The  three  windows  used  by  the  AFGWC  Boundary-layer  Model  are  also  shown  in 
Fig.  3.23  (inner  dashed  lines)  and  summarized  in  Table  3.6.  A  29x27  half-mesh 
window  for  the  U.S.  and  29x35  half-mesh  windows  for  Europe  and  Asia  art* 
defined. 


3. 1.3.7.  MOVABLE  WINDOW  GRIDS 

One  set  of  movable  window  grids  is  used  at  AFGWC  by  the  II igh-reso lut ion 
Cloud  Prognosis  Model  (HRCP).  These  grids  are  movable  in  the  sense  that  some 
flexibility  exists  in  specifying  their  location.  Each  window  grid  is  one  of 
the  3DNEPH  boxes,  discussed  in  Section  3. 1.3.3,  and  contains  an  array  of  64x69 
eighth-mesh  grid  points.  Fig.  3.17  provides  an  example  of  one  such  box  (box 
number  44). 

The  particular  3DNEPH  boxes  are  selected  according  to  areas  of  interest 
and  may  be  changed  each  time  HRCP  is  used.  Up  to  12  boxes  may  be  specified, 
although  typically  about  8  are  selected.  Forecasts  for  the  boxes  arc  made 
independently  of  each  other  by  HRCP. 
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Half-mesh  Analysis  Model  grids  (solid  lines)  ai 
Boundary-layer  Model  grids  (inner  dashed  lines] 
domain  of  the  Northern  Hemispheric  Reference  Gi 
outer  border)  and  the  Northern  Hemispheric  Octi 
(outer  dashed  lines)  are  indicated.  The  indic< 
(IW,JW)  and  (I2.J2)  designate  the  coordinates  1 
the  whole-mesh  and  half-mesh  versions,  respect: 
respect  to  the  Octagons  and  the  Reference  Gridi 
whole-mesh  grid  spacing  is  displayed  in  the  up] 


corner. 


3.2  MERCATOR  GRIDS 


Mercator  grids  are  used  at  AFGWC  for  meteorological  applications  in 
tropical  regions.  There  are  two  such  grids.  One  we  designate  as  the 
Conventional  Tropical  Grid  (Fig  3.24),  which  is  the  grid  used  for  the 
conventional  meteorological  elements,  such  as  height,  wind,  and  temperature. 
The  second  grid  is  the  SGDB  Tropical  Grid  (Fig.  3.25),  which  is  used  for 
processing  satellite  imagery.  Both  grids  are  on  a  Mercator  projection  true  at 
22.5°N  and  22.5°S.  The  Conventional  Tropical  Grid  is  73x19  with  a  grid 
length  of  5°  longitude,  whereas  the  SGDB  Whole-mesh  Tropical  Grid  is  97x28 
with  a  grid  length  of  about  3.71°  longitude.  The  Greenwich  Meridian  is 
represented  by  columns  1  and  73  in  the  former,  but  only  by  column  1  in  the 
latter.  If  the  SGDB  Whole-mesh  Tropical  Grid  had  a  98th  column,  this  column 
would  be  coincident  with  the  Greenwich  Meridian. 

The  longitudinal  (east-west)  spacing  AA  between  grid  points  is  given  by 

360°/T  ,  in  degrri , 

o  /t  0-35) 

<?7t/T^  ,  in  radian:., 

where  Ijq  is  the  number  of  nonduplicated  points  in  the  east-west  direction. 

For  the  Conventional  Tropical  Grid  1^  =  72,  so  that  AA  =  5°  =  0.087266 
radians  longitude.  Because  IN  =  97  for  the  SGDB  Tropical  Grid,  AA  is  about 
3.7113°  or  0.064775  radians  longitude. 

The  distance  d^  between  grid  points  on  the  image  plane  for  a  Mercator 
projection  is  given  by 

d.  =  a  (cos  4.  )  AA  ,  (3.36) 

where  a  is  the  earth's  radius,  <f  ^  is  a  standard  latitude  (2  2 . 5°N  at  AFGWC), 
and  AA  is  the  longitudinal  spacing  in  radians  calculated  from  eq.  (3.35).  As 
mentioned  above  for  the  Conventional  Tropical  Grid,  the  spacing  is  5°  = 
0.087266  radians  longitude,  so  that  according  to  eq.(3.36),  for  =  22.5°N, 
we  obtain  d^  =  513.67...  km.  For  the  SGDB  Whole-mesh  Tropical  Grid,  AA  is 
about  3.7113°  =  0.064775  radians,  so  tiiat  d^  =  381.28...  km. 

The  distance  de  between  grid  points  on  the  earth's  surface  can  be 
derived  from  d^  by  applying  the  definition  of  image  scale  0  (Section  2.2)  to 
eq.  (3.36): 

d. 

d  =  —  =  a  (cos  1(1  )  AA  ,  (3.37) 

e  0 

where  <j  is  specified  by  eq.  (2.9)  using  standard  latitude  <l>j  of  22.5°N,  41  is 
latitude,  and  AA  is  the  longitudinal  specing  in  radians  (see  preceding 
paragraph).  The  variation  of  distance  de  between  grid  points  is  shown  ns  a 
function  of  latitude  in  Fig.  3.26.  Because  the  grid  length  de  for  the 
whole-mesh  polar  stereographic  grids  and  for  the  SGDB  Whole-mesh  Tropical  Grid 
are  approximately  equal,  the  la'. ter  is  designated  the  whole-mesh  grid  length 
for  Mercator  projections. 
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Y  =  -  (J  -  J£)  ,  1  i  J  i  JM  -  (3. AO) 

In  these  equations  K  is  any  integer,  although  usually  0  or  -1.  Jg  is  the  J 
coordinate  of  the  Equator.  Ij^  is  the  total  number  of  points  in  the 
east-west  direction.  is  the  number  of  nonduplicated  points  in  the 

east-west  direction,  with  point  1=1  being  coincident  with  point  1^+1. 

JM  is  the  number  of  points  in  the  north-south  direction.  Finally,  G 
designates  the  grid,  with  G  =  -1  for  the  Conventional  Tropical  Grid  and  G  =  +1 
for  the  SGDB  Whole-mesh  Tropical  Grid.  Table  3.7  contains  a  summary  of  the 
values  of  these  parameters  for  the  two  Mercator  grids. 


Because  the  grid  length  de  is  a  function  of  latitude,  the  distance  de' 
between  grid  points  on  the  earth's  surface  in  the  north-south  direction  is  not 
quite  equal  to  de.  Instead,  we  have  that 


d'  =  a  ;,{  ,  (3.38) 

e 

where  A4>  is  the  latitudinal  spacing  (in  radians)  between  adjacent  grid 
points.  As  stated  above,  the  distance  dg  between  adjacent  grid  points  along 
a  given  latitude  corresponds  to  5°  longitude  for  the  Conventional  Tropical 
Grid.  By  equating  de  and  de*  in  eqs.  (3.37)  and  (3.38),  we  can  show  that 
distance  de 1  between  adjacent  grid  points  along  a  meridian  corresponds  to 
about  3°  latitude  at  the  Equator  and  about  3.8°  latitude  near  the  northern 
and  southern  boundaries.  Similarly,  for  the  SGDB  Whole-mesh  Tropical  Grid, 
for  which  de  corresponds  to  about  3.7°  longitude,  de‘  corresponds  to 
about  3.7°  latitude  at  the  Equator  and  about  2.6°  latitude  near  the 
northern  and  southern  boundaries. 


The  (I,J)  indexing  conventions  for  the  two  grids  are  shown  in  Figs. 3. 24 
and  3.25.  The  convention  for  the  Conventional  Tropical  Grid  is  right-handed; 
that  is,  I  increases  from  right  to  left  (east  to  west)  and  J  increases  from 
top  to  bottom  (north  to  south).  The  convention  for  the  SGDB  Whole-mesh 
Tropical  Grid,  however,  is  left-handed;  I  increases  from  west  Lo  east  and  J 
increases  from  north  to  south. 


Grid-point  Cartesian  coordinate  systems  (X,Y)  are  defined  for  botli 
tropical  Mercator  grids.  The  system  is  right-handed  for  both  (Figs. 3.24  and 
3.25).  The  origin  for  both  systems  is  located  at  the  intersection  of  the 
Equator  and  the  Greenwich  Meridian.  Equations  for  (X,Y)  in  multiples  of 
distance  d^  between  grid  points  on  the  image  plane  are 


Table  3.7.  Parameters  for  the  Mercator  Grids 


PARAMETER 

SYMBOL 

CONVENTIONAL 

TROPICAL  GRID 

SGDB  WHOLE-MESH 
TROPICAL  GRID 

Standard 

latitudes 

*1'*2 

22 . 5°N  ,  22 . 5°S 

22.5°N,  22 . 5°S 

Reference 

longitude 

* 

0 

0°E 

0°E 

Longitudinal 
grid  spacing 

A  A 

5°  longitude 

3.7113°  longitude 

Distance  between 
grid  points  on 
the  image  plane 

di 

513.67  km 

381.28  km 

(I,J)  axes 

right-handed 

left-handed 

(X,Y)  axes 

right-handed 

right-handed 

Total  number  of 
points  in  the 
east-west 

direction 

IM 

73 

97 

Number  of 

nonduplicated 
points  iu  the 
east-west 
di rection 

XN 

72 

97 

Number  of 

points  in  the 

north-south 

direction 

JM 

19 

28 

Range  of  I  index 

1  i  I  1  IM 

1  1  I  1  IN+1 

1  H  llM 

1  11  llN 

Range  of  J  index 

1  1  J  1  JM 

l  1  J  1  Jm 

J  coordinate  of 
the  Equator 

JE 

10 

15 

Grid  designator 

G 

-1 

+  1 

Image- plane  Cartesian  coordinates  (x,y)  can  be  calculated  from  (X,Y) 
coordinates  by  multiplying  each  (X,Y)  value  by  d^  (the  image-plane  distance 
between  grid  points): 


x  =  d .  X 

i 


(3.41) 


Y  =  d.  Y 


(3.42) 


Using  eqs.  (2.10),  (2.11),  (3.36),  and  ( 3 . 39  )- ( 3 . 42  )  ,  we  can  derive  the 
latitude  <p  and  longitude  A  of  any  grid  point  from  its  (I,J)  coordinates: 


i(>  =  ;  arc  tan  -{exp 


•(  AA  )  (,l  -.1 ) 


(3.43) 


0(1-1 ) ( AA )  +  \ 

o 

G  ( 1-1  -1  )  ( AA  )  +  A 
N  o 


(3.44) 


where  the  various  parameters  are  provided  in  Table  3.7.  In  eq.  (3.43)  note 
that  AA  and  the  functional  value  of  the  arctangent  should  be  in  degrees. 

The  (I,J)  coordinates  of  any  point  within  the  latitudinal  bounds  of  the 
tropical  grids  may  be  calculated  from  latitude  and  longitude  (eqs.  (3.43)  and 
(3.44))  according  to 


SGDB  WhO-Lo-munli  Tropical  Hr  id  , 


(3.45) 


1  +1  - 
N 


Conventional  Tropical  Grid 


J  =  JE  ' 


u  ( AA  ) 


£n[tun  hi  <t>+90° )  1 


(3.46) 


In  eq.  (3.45)  A  must  have  a  range  of  \Q  A  c  Ao+360°  for  the  SGDB 
Whole-mesh  Tropical  Grid  and  A0  <_  A  <_  AQ  +360°  for  the  Conventional  Tropical 
Grid. 


3.2.1  Conventional  Tropical  Grid 

The  Conventional  Tropical  Grid  is  formed  by  superimposing  a  73x19  array  of 
points  on  a  Mercator  projection  such  that  row  10  (J  =  10)  coincides  with  the 
Equator  and  columns  1  (I  °  1)  and  73  (I  “  73)  both  coincide  with  the  Greenwich 
Meridian.  As  shown  in  Fig.  3.24,  the  grid  extends  completely  around  the  world 
longitudinally  and  includes  that  portion  of  the  globe  between  approximately 
41°N  and  41°S.  This  grid  is  one  of  the  few  at  AFGWC  to  use  a  right-handed 
(I,J)  indexing  convention.  Equations  governing  this  grid  are  presented  in  the 
preceding  section,  with  the  relevant  grid  parameters  specified  in  Table  3.7. 

The  grid-point  Cartesian  coordinate  system  (X,Y)  provides  the  frame  of 
reference  for  orientation  of  wind  components:  the  u  component  is  positive  in 
the  direction  of  decreasing  X  and  the  v  component  is  positive  in  the  direction 
o f  decreasing  Y  (Fig.  3.24). 

3.2.2.  SGDB  Tropical  Grid  System 

The  SGDB  Whole-mesh  Tropical  Grid  is  formed  by  superimposing  a  97x28  array 
of  points  on  a  Mercator  projection  such  that  row  15  (J  =  15)  coincides  with 
the  Equator  and  column  1  coincides  with  the  Greenwich  Meridian.  If  a  98th 
column  were  included,  it  would  also  coincide  with  the  Greenwich  Meridian.  As 
shown  in  Fig.  3.25  the  grid  extends  completely  around  the  world  longitudinally 
and  from  about  46°N  to  46°S  lati tudinal ly .  Equations  governing  this  grid 
are  presented  in  Section  3.2,  with  the  relevant  grid  parameters  specified  in 
Table  3.7. 

As  was  similarly  true  for  SGDB  grids  on  polar  stereographic  projections, 
the  SGDB  Whole-mesh  Tropical  Grid  is  partitioned  into  boxeB  to  facilitate 
processing  of  satellite  data.  Each  SGDB  box  contains  one  grid  point 
(*w>Jw)  of  the  SGDB  Whole-mesh  Tropical  Grid  and  consists  of  a  64x64  array 
of  sixty-fourth-mesh  ooints.  There  are  97x28  =  2716  boxes  containing  a  total 
of  11,124,736  grid  points.  Each  point  UW,JW)  coincides  with  the 
sixty-fourth-mesh  point  in  the  upper-left  corner  of  the  SGDB  box.  SGDB  boxes 
are  numbered  according  to 


B  =  (J  -1)  +  -1) 

W  W 


where  (IW,JW)  are  the  indices  for  the  SGDB  Whole-mesh  Tropical  Grid  of  the 
point  at  the  upper-left  corner  of  the  box  and  Bs  is  the  SGDB  box  number. 
Note  that  0  Bs  2715. 


1  ■ 

1  <  ,l  •  i'8  , 

—  w  — 


(3.47) 
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Satellite  data  processing  software  (SYNAPS)  designates  the  origin  of  the 
SGDB  Whole-mesh  Tropical  Grid  as  the  point  (0,0),  as  opposed  to  (1,1),  the 
convention  of  this  Technical  Note.  Consequently, 

I '  =  1  -  1  ,  (3.48) 

w  w 


,]'=.!  -  1 
W  w 


where  (Iw' ,  Jw' )  are  the  whole-mesh  indices  as  used  in  SYNAPS. 


(3.49) 


3.3  LATITUDE-LONGITUDE  GRID 


The  Global  Applications  Data  Base  (GADB)  grid  is  the  only 
latitude-longitude  grid  currently  in  use  at  AFGWC.  It  consists  of  a  73x52 
array  of  points  extending  around  the  world  in  a  band  between  87.5°N  and 
87.5°S.  This  grid  is  formed  by  allowing  each  point  to  represent  a  specific 
latitude  and  longitude  (Fig.  3.27).  For  example,  grid  point  (1,1),  in  the 
upper-left  corner  of  the  grid,  represents  2 . 5°E/87 . 5°N .  Grid  spacing  is 
5°  longitude  in  the  east-west  direction;  this  longitudinal  spacing  is  the 
same  as  the  Conventional  Tropical  Grid.  In  the  north-south  direction,  spacing 
is  2.5°  latitude  in  the  region  17.5°N  to  57.5°N  and  17.5°S  to 
57,5°S,  and  5°  latitude  elsewhere. 

For  this  grid  the  (l,J)  indexing  convention  is  left-handed;  that  is,  the  I 
index  increases  from  left  to  right  (west  to  east)  as  J  increases  from  top  to 
bottom  (north  to  south).  Columns  I  B  1  and  1  =  73  represent  the  same 
longitude,  2.5°E,  with  the  Greenwich  Meridian  being  centered  between  columns 
I  =  72  and  I  =  73.  The  Equator  lies  midway  between  rows  J  =  26  and  J  =  27. 


SECTION  4.  APPENDIX  A— THE  REGIONS  OF  THE  AFGWC  REGIONS  DATA  BASE 


World-wide  observations  from  land  surface  stations  and  radiosondes,  for 
example,  are  often  located  geographically  by  the  WMO  according  to  WMO 
Regions .  Figs.  A. 1  and  A. 2  indicate  the  WMO  Regions  for  the  Northern  and 
Southern  Hemispheres,  respectively.  Because  these  regions  are  occasionally 
modified,  an  up-to-date,  authoritative  reference  should  be  consulted  when 
knowledge  of  exact  regions  is  important. 

At  AFGWC,  these  observations  are  not  stored  at  grid  points  but  rather  are 
stored  individually  by  geographic  regions  in  the  so-called  Regions  Data  Base. 
Fig.  A. 3  shows  the  domains  of  the  65  Northern  Hemispheric  regions  and  Fig.  A. 4 
the  domains  for  the  11  Southern  Hemispheric  regions. 


gure  A.  1.  World  Meteorological  Organization  Regions  for  lu  ating 

meteorological  observations  in  the  Northern  Hemisphere. 
The  domain  of  the  Northern  Hemispheric  Reference  grids  is 
indicated  by  the  solid  outer  border. 
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Woria  Meteorological  urganization  Regions  ror  locating 
meteorological  observations  in  the  Southern  Hemisphere. 
The  domain  of  the  Southern  Hemispheric  Reference  grids  is 
indicated  by  the  solid  outer  border. 


AFGWC  Southern  Hemispheric  Regions  for  storing  observations 
in  the  Regions  Data  Base.  There  are  11  Southern 
Hemispheric  Regions,  numbered  from  bb  to  7b.  The  domain  of 
the  Southern  Hemispheric  Reference  Grids  is  indicated  by 
the  solid  outer  border. 


SECTION  5.  APPENDIX  B  -GLOSSARY 


This  is  an  alphabetical  listing  of  definitions  and  acronyms  used  in  this 
Technical  Note-  A  List  of  Symbols  is  presented  elsewhere.  Location  of 
additional  information  or  primary  usage  is  indicated,  when  appropriate. 


AFGWC 


-  Air  Force  Global  Weather  Central. 


AGL 

Analysis 


Asian  boundary- 
layer  window 

Asian  window 


AWSPE 


Bedient 

Boundary-layer  Model 


Class 


Coincidence  class 


Conformal  variety 


Conical  variety 


Constant  of  the 
cone  (n) 


-  Above  ground  level. 

-  Process  in  which  data  are  interpolated  from  the 
observational  points  to  the  grid  points  (Section  3). 

-  Areal  subset  of  the  Asian  window;  used  by  the  Asian 
Boundary-layer  Model  (Section  3. 1.3. 6. 2). 

-  Stationary  window  covering  eastern  Asia;  used  by  the 
Asian  Half-mesh  Analysis  Models  (Section  3. 1.3. 6.1). 

-  Air  Weather  Service  Primitive  Equation  Model; 
primary  whole-mesh  forecast  model  for  the  Northern 
Hemisphere;  provides  forecasts  from  the  earth's 
surface  to  100  mb  (Section  3. 1.3. 2). 

-  Same  as  whole-mesh  grid  length  (Section  3.1.1). 

-  Model  that  provides  half-mesh  forecast',  for  eight 
levels  from  the  earth's  surface  to  1600  m  AGL;  used 
with  the  Asian,  European,  and  1J.S.  boundary-layer 
windows  (Section  3. 1.3. 6. 2). 

-  One  of  five  different  categories  by  which 
projections  can  be  grouped.  The  classes  are  not 
mutually  exclusive  (Section  2.3). 

-  Indicates  the  type  of  contact  between  the  earth  and 
the  geometric  figure  associated  with  a  projection 
(Section  2.3). 

-  A  variety  of  the  property  class  in  which  the  angles 
between  intersecting  curves  and,  therefore,  the 
shapes  of  small  features  on  the  earth's  surface  are 
maintained.  A  differentially  small  distance  on  the 
projection  will  correspond  to  a  distance  on  the 
earth's  sjrface  independent  of  direction  from  a 
point  (Section  2.3). 

-  /..  variety  of  the  nature  class  in  which  a  cone  is  the 
geometric  figure  associated  with  the  projection 
(Section  2.3). 

-  A  scaling  factor  used  in  Lambert  conformal 
projections  (see  Section  2.6.2,  eq.  (2.13)). 
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Conventional  Tropical 
Grid 

Conventional  variety 

Cylindrical  variety 

De ve loping 

Domain 

Eighth-mesh  grid 
Equidistant  variety 

Equivalent  variety 

European  boundary- 
layer  window 

European  window 

Finer-mesh  grid 

GADB 

Generation  class 

Geometric  variety 

Geometric  figure 


-  A  tropical  grid  on  a  Mercator  projection  used  for 
the  conventional  meteorological  elements,  such  as 
height,  wind,  and  temperature  (Section  3.2). 

-  A  variety  of  the  generation  class  in  which  the 
projection  is  derived  from  the  earth's  surface 
mathematically,  without  using  rays  (Section  2.3). 

-  A  variety  of  the  nature  class  in  which  a  cylinder  is 
the  geometric  figure  associated  with  the  projection 
(Section  2.3). 

-  The  process  of  cutting  and  unfolding  a 
three-dimensional  surface  to  produce  a  planar 
surface  (Section  2.2). 

-  The  horizontal  area  represented  by,  for  example,  a 
map  or  grid. 

-  Grid  on  which  the  grid  length  is  one-eighth  the 
whole-mesh  value. 

-  A  variety  of  the  property  class  in  which  distances 
between  points  on  the  earth  are  correctly 
represented  on  the  projection  (Section  2.3). 

-  A  variety  of  the  property  class  in  which  areas  of 
features  on  the  earth  are  correctly  represented  on 
the  projection  (Section  2.3). 

-  Areal  subset  of  the  European  window;  used  by  the 
European  Boundary-layer  Model  (Section  3. 1.3. 6. 2). 

-  Stationary  window  covering  Europe  and  western  Asia; 
used  by  the  North  American  Half-mesh  Analysis  Models 
( Section  3 . 1. 3 .6 . 1) . 

-  Grid  with  a  grid  length  shorter  than  the  whole-mesh 
grid  length  (Section  3). 

-  Global  Applications  Data  Base;  data  base  containing 
global  analyses  and  forecasts  of  various 
meteorological  elements  used  primarily  in  computer 
flight  planning  (Section  3.3). 

-  Indicates  the  method  used  to  construct  the  image 
surface  in  generating  a  projection  (Section  2.3). 

-  A  variety  of  the  generation  class  in  which  the 
surface  of  the  earth  is  simply  projected  using  rays 
onto  a  geometric  figure  in  generating  a  projection 
(Section  2.3). 

-  Two-  or  three-dimensional  surface  associated  with  a 
projection  of  the  earth's  surface  (Section  2.2). 
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Great  circle 


Grid 

Grid  increment 

Grid  interval 

Grid  length 

Grid-point  Cartesian 
coordinates  (X,Y) 

Grid  unit 
Half-mesh  grid 

Half-mesh  Analysis 
Model 

HRCP 

HUFANL 

IJ  indexing  conventioi 
Image  plane 

Image-plane  Cartesian 
coordinates  (x,y) 

Image-plane  polar 
coordinates  (r,0) 

Image  scale  (a) 

Image  surface 


-  A  line  of  shortest  length  between  two  points  along 
the  earth's  surface  (Section  2.6.1). 

-  Array  of  points  used  to  represent  locations  on  the 
surface  of  the  earth  (Section  3). 

-  Same  as  grid  length. 

-  Same  as  grid  length. 

-  Distance  between  adjacent  points  on  a  grid. 

-  A  Cartesian  coordinate  system  for  a  grid  in  which 
locations  are  specified  by  multiples  of  the  grid 
length  for  that  particular  grid  on  the  image  plane. 

-  Same  as  grid  length. 

-  Grid  on  which  the  grid  length  is  one-hal'"  the 
whole-mesh  value. 

-  Numerical  weather  analysis  model  which  uses  a 
half-mesh  grid. 

-  High-resolution  Cloud  Prognosis  Model;  cloud 
forecast  model  with  fifteen  levels  in  the  vertical 
and  eighth-mesh  horizontal  resolution  (Section 

3. 1.3.7). 

-  Hough  Global  Spectral  Analysis  Model;  provides 
whole-mesh  global  analyses  for  the  mandatory 
pressure  levels  from  the  earth's  surface  to  100  mb. 

-  System  for  uniquely  identifying  individual  grid 
points. 

-  The  two-dimensional  surface  to  which  the  surface  of 
the  earth  is  transformed  by  projection,  by  pure 
mathematical  equations,  or  by  a  combination  of  both 
a  planar  image  surface  or  an  image  surface  that  has 
been  developed  into  a  plane  (Section  2.2). 

-  A  Cartesian  coordinate  system  that  is  on  the  image 
plane  and  is  independent  of  grid  length. 


-  \  polar  coordinate  system  on  the  image  plane. 


-  Ratio,  at  a  point,  of  differential  distance  on  an 
image  surface  or  image  plane  to  differential 
distance  on  the  surface  of  the  earth  (Section  2.2). 

-  The  surface  (not  necessarily  two-dimensional)  to 
which  the  surface  of  the  earth  is  transformed  by 
projection,  by  pure  mathematical  equations,  or  by  a 
combination  of  both  (Section  2.2). 
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Mathematical  function  whose  functional  value  is  the 
integer  portion  of  the  division  of  the  first 
argument  by  the  second  argument;  for  example, 

INT( 10, 3 )  =  3  (Section  3. 1.3.3). 


Lambert  conformal 
projection 


Longitudinal  grid 
spacing 

Loxodrome 


Map  factor  (m) 


Map  scale  (y) 


Mercator  projection 


Meridian 


Mesh  factor  (M) 


Mesh  size 


Movable  window 
grid 

Nature  class 


One  type  of  conformal  projection  using  a  cone  as  the 
image  surface;  common  in  limited-area  raidlatitudinal 
applications  (Section  2.6). 

Spacing  of  grid  points  in  the  east-west  direction. 


Line  or  path  that  maintains  fixed  compass  direction; 
shown  on  a  map  as  a  line  crossing  all  meridians  at 
the  same  angle;  also  called  rhumb  line  (Section 
2.5.1). 

Two-dimensional  horizontal  representation  of  the 
surface  of  the  earth  with  size  convenient  for 
display,  analysis,  and  interpretation  of  information 
(Section  2.2). 

Ratio,  at  a  point,  of  differential  distance  on  a  map 
to  differential  distance  on  the  earth;  the  product 
of  image  scale  o  and  map  scale  V  (Section  2.2). 

Ratio  of  distance  on  a  map  to  distance  on  the  image 
surface  or  image  plane;  describes  mathematically  the 
uniform  reduction  from  image  surface  or  image  plane 
to  map  surface  (Section  2.2). 

One  type  of  conformal  projection  using  the  cylinder 
as  the  image  surface;  common  in  tropical 
applications  (Section  2.5). 

Line  of  constant  longitude,  extending  from  the  North 
Pole  to  the  South  Pole. 

Ratio  of  the  grid  length  to  the  whole-mesh  grid 
length  on  the  image  plane  (Section  3.1.1). 

Indicator  of  the  relative  grid  length  for  a 
particular  grid  (for  example,  whole  mesh,  half  mesh, 
eighth  mesh)  (Section  3. 1.2. 4). 

Mathematical  function  in  which  the  functional  value 
is  the  remainder  of  the  division  of  the  first 
argument  by  the  second  argument;  for  example, 

M0D( 10, 3)  =  1  (Section  3. 1.3. 3). 

Grid  with  a  limited-area  domain  whose  location 
can  be  varied  (Section  3. 1.3. 7). 

Indicates  the  geometric  figure  associated  with  a 
projection  (Section  2.3). 
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NH  -  Northern  Hemisphere  (Table  3.5). 

Normal  variety  -  A  variety  of  the  position  class  in  which  the  axis  of 

symmetry  of  the  geometric  figure  associated  with  the 
projection  coincides  with  the  rotational  axis  of  the 
earth  (Section  2.3). 


North  American  window  -  Stationary  window  covering  most  of  North  America; 

used  by  the  North  American  Half-mesh  Analysis  Models 
(Section  3. 1.3.6. 1). 

Oblique  variety  -  A  variety  of  the  position  class  in  which  the  axis  of 

symmetry  of  the  geometric  figure  associated  with  the 
projection  is  neither  perpendicular  to  nor 
coincident  with  the  rotational  axis  of  the  earth 
(Section  2.3). 


Octagon 
Octagon  grid 


Orthomorphic 

variety 


-  Same  as  Octagon  Grid. 

-  An  octagon-shaped  grid;  an  important  numerical 
weather  analyses  and  prediction  grid  for  both  the 
Northern  and  Southern  Hemispheres  (Section  3. 1.3.1). 

-  Same  as  conformal  variety. 


Parallel 


-  Line  of  constant  latitude,  encircling  the  earth  in 
the  east-west  direction. 


Perspective  variety  -  Same  as  geometric  variety. 


Planar 


-  Pertaining  to  a  plane;  two-dimensional. 


Planar  variety 


-  A  variety  of  the  nature  class  in  which  a  plane  is 
the  geometric  figure  associated  with  the  projection 
(Section  2.3). 


Polar  stereographic 
projection 


-  One  type  of  conformal  projection  using  a  plane 
as  the  image  surface;  common  in  Northern  and 
Southern  Hemispheric  applications  (Section  2.4). 


Polysuperf icial  -  A  variety  of  the  coincidence  class  in  which  a 

variety  succession  of  different  geometric  figures  is 

associated  with  the  projection  (Section  2.3). 


Position  class 


Projection 


-  Indicates  the  alignment  of  the  axis  of  the  earth 
with  respect  to  the  geometric  figure  associated  with 
a  projection  (Section  2.3). 

-  The  representation  of  the  surface  of  the  earth  on  an 
image  surface  or  image  plane.  In  some  cases  this 
representation  can  be  generated  geometrically  using 
rays  from  a  point  on  or  within  the  earth's  surface 
to  project  the  earth's  surface  onto  a  geometric 
figure  (Section  2.2). 
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Property  class 


-  Indicates  the  conservation  of  distance,  area,  shape 
or  angle  in  generating  a  projection  (Section  2.3). 


Quarter-raesh  grid 
Ray 

Reference  grid 

Reference  longitude 

Regions  Data  Base 

Rhumb  line 
Secant  variety 

Semi-geometric 

SGDB 

SGDB  box 

SGDB  Box  IJ 
convention 

SGDB  Box  IJ 
indices 

SGDB  Tropical 
Grid 

SH 

SIXLVL 

Sixty- fourth-mesh 
grid 


-  Grid  on  which  the  grid  length  is  one-foui, th  the 
whole-mesh  value. 

-  Straight  line  used  to  project  the  surface  of  the 
earth  onto  the  geometric  figure. 

-  The  areal  superse1-  grid  for  all  grids  of  a  given 
mesh  size  and  hemisphere. 

-  Longitude  used  in  specifying  the  coordinate  axes  for 
projections,  maps,  and  grids. 

-  AFGWC  data  base  in  which  observations  are  stored 
individually  by  geographic  regions  (Section  A). 

-  Same  as  loxodrome. 

-  A  variety  of  the  coincidence  class  in  which  the 
geometric  figure  associated  with  the  projection 
intersects  the  earth's  interior  (Section  2.3). 

-  A  variety  of  the  generation  class  in  which  the 
projection  is  derived  from  the  earth's  surface  using 
a  combination  of  projection  by  rays  and  mathematics 
without  consideration  of  rays  (Section  2.3). 

-  Satellite  Global  Data  Base;  data  base  that  contains 
global  high-resolution  (about  3  nm  or  6  km)  radiance 
measurements  (from  satellite)  in  the  visual  and 
infrared  ranges  (Sections  3. 1.3.5,  3.2). 

-  Areal  subset  of  the  SGDB  grid;  each  box  contains 
64x64  =  4096  sixty-fourth-mesh  points  (Sections 
3. 1.3. 5,  3.2.2). 

-  Convention  for  uniquely  specifying  each  of  the 
64x64  grid  points  in  a  SGDB  box  (Section  3. 1.3.5). 

-  The  indices  of  the  SGDB  Box  IJ  convention 
(Section  3 . 1 . 3 . 5 ) . 

-  A  tropical  grid  on  a  Mercator  projection  used 
for  processing  satellite  imagery  (Section  2.3). 

-  Southern  Hemisphere  (Table  3.5.). 

-  Six-level  Baroclinic  Prediction  Model;  provides 
operational  forecasts  for  the  Southern  Hemisphere 
and  updated  forecasts  for  the  Northern  Hemisphere 
(Section  3. 1.3. 2). 

-  Grid  on  which  the  grid  length  is  one-sixty-fourth  the 
whole-mesh  value. 


Standard  latitude(s) 


Stationary  window 
grid 

Super  Grid 
SYNAPS 


Tangent  variety 

Transverse  variety 

TRONEW 

True  latitude 

U.S.  boundary-layer 
window 

U.S.  window 
Variety 

Whole-mesh  grid 

Whole-mesh  grid 
length 

Whole-mesh  grid 
spacing 


-  Latitude(s)  at  which  the  image  surface  represents 
true  earth  distance;  that  is,  latitude(s)  at  which 
the  image  scale  a  is  unity. 

-  Grid  with  a  limited-area  domain  whose  location 
cannot  lie  varied  (Section  J. 1.3.6). 

-  Same  as  Reference  Grid. 

-  Program  that  processes  video  and  infrared  satellite 
data  and  stores  them  into  the  SGDB;  the  processing 
interface  between  the  satellite  data  fonnaters  and 
the  SGDB.  The  name  is  derived  from  the  word 
synapse,  which  is  the  point  of  contact  along  which 
impulses  are  transmitted  between  two  nerve  cells 
(Sections  3. 1.3. 5,  3.2.2). 

-  A  variety  of  the  coincidence  class  in  which  the 
geometric  figure  associated  with  the  projection 
intersects  only  the  earth's  surface  and  not  the 
earth's  interior  (Section  2.3). 

-  A  variety  of  the  position  class  in  whicli  the  axis  of 
symmetry  of  the  geometric  figure  associated  with  the 
projection  is  perpendicular  to  the  rotational  axis 
of  the  earth  (Section  2.3). 

-  Tropical  Cloud  Forecast  Model  (Section  3. 1.3.4). 

-  Same  as  standard  latitude. 

-  Areal  subset  of  the  North  American  window; 
used  by  the  U.S.  Boundary-layer  Model  (Section 
3.  1.3. 6. 2). 

-  Same  as  North  American  window. 

-  One  of  the  mutually  exclusive  subsets  of  each  of  the 
five  classes  used  in  categorizing  projections 
(Section  2.3). 

-  Grid  whose  grid  length  serves  as  a  basis  tor  /ill 
finer-meali  grid  lengths;  grid  whose  grid  length  in 
the  whole-mesh  grid  length  (Sections  l.l.l,  3.2). 

-  For  a  polar  stercogruphic  projection,  u  grid 
length  of  exactly  381  km  at  60°N  (or  °s)  on  the 
earth's  surface;  for  a  Mercator  projection,  a  grid 
length  of  approximately  381  km  at  22.3°N  and  °S 

on  the  earth's  surface  (Sections  3.1.1,  3.2). 

-  Same  as  whole-mesh  grid  length. 


WMO 


-  World  Meteorological  Organization. 


WMO  Region 


3DNEPH 


3DNEPH  box 


3DNEPH  Box  IJ 
convention 

3DNEPH  Box  IJ 
indices 


-  A  geographic  region  by  which  the  location  of 
observations  may  be  specified  according  to 
convention  established  by  the  WMO  (Section  A) . 

-  Three-dimensional  Nephanalysis  Model;  provides 
eighth-mesh  cloud  analyses  for  both  the  Northern  and 
Southern  Hemispheres  (Section  3. 1.3. 3). 

-  One  of  64  areal  subsets  of  the  3DNEPH  grid;  each  box 
contains  64x64  n  4096  eighth-mesh  points  (Section 

3. 1.3. 3). 

-  Convention  for  uniquely  specifying  each  of  the 
64x64  grid  points  in  a  3DNEPH  box  (Section  3. 1.3.3). 

-  The  indices  of  the  3DNEPH  Box  IJ  convention 
(Section  3. 1 . 3 .3 ) . 
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